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Pref  ace 


The  purpose  of  this  study  was  to  verify  the 
mathematical  correctness  and  to  determine  some  of  the 
limitations  of  the  monostati c-bi static  radar  cross  section 
(RTS)  relationship  developed  by  Robert.  E.  Kell. 

Pot ent i a  1 1 y , t he  Radar  Target  Scattering  Facility,  Holloman 
APR,  NM ,  could  apply  this  re 1  a t i onsh i p  to  reduce  testing 
time  and  costs. 

-\  detailed  mathematical  formulation  is  provided.  It 
begins  with  the  monostat i c-bi s tat i c  equivalence  theorem 
( MRFT ) .  The  MBET  is  a  physical  optics  approximation 
relating  monostatic  and  bistatic  RCS .  The  derivation  of 
Kell’s  relationship  is  built  upon  the  MBET  and  is 
mathematically  sound  within  the  prescribed  bounds. 

Kell's  relationship  was  tested  using  computer 
models.  The  results  of  the  testing  show  Kell’s  method  of 
determining  bistatio  RCS  from  monostat. ic  RCS  measurements 
has  some  merit..  Measurements  are  needed  to  establish  the 
limits  to  which  Kell’s  method  may  be  applied.  Also,  the 
complexity  of  the  targets  mu9t  be  increased  in  order  to  gage 
the  accuracy  of  Kell’s  method  under  a  variety  of  conditions. 

In  performing  this  study  1  have  received  help  from 
others.  Of  course,  I  am  indebted  to  my  faculty  advisor, 
Major  Harry  Barksdale,  whose  insight  and  guidance  kept  me 
moving  in  the  right  direction.  I  would  also  like  to  thank 
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Captain  Ri  1 1  Oet.ting  who  never  failed  to  respond  positively 
to  my  panicked  cries  for  help  and  Captain  Dennis  Tackett  who 
devoted  a  great  deal  of  time  and  effort  in  the  actual 
production  of  this  thesis.  And,  finally,  a  special 
thanks  to  my  wife,  DeV»bie,  who  endured  it.  all  with  a  smile. 
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Robert,  Kell  developed  a  relationship  between  monostatic 
radar  cross  section  (RCS)  and  bistatic  RCS  whereby  bistat.io 
RCS  can  be  predicted  from  monost.atic  RCS  measurements  under 
certain  conditions.  This  study  found  Kell's  relationship  to 
be  mathematically  sound  given  certain  assumptions.  Kell's 
relationship  was  then  tested  by  comparing  computer  generated 
bistatic  and  monostntic  cross  sections  for  simple  shapes. 
Four  parameters  were  varied  during  testing  in  order  to 
discern  possible  limitations  of  Kell's  method:  bistatic 
angle,  angle  of  incidence,  electrical  siz  of  the  target, 
and  continuity.  Results  of  the  testing  show  Kell's  method 
has  some  merit.  The  difference  between  the  bistatic  RCS  and 
its  related  monost.atic  RCS  for  electrically  large  spheres  is 
less  than  1-dB  up  to  bistatic  angles  of  80*^  For 
electrically  large  flat  and  singly  curved  surfaces  the 
monost.atic  and  bistatic  cross  sections  were  within  3-dB  for 
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angles  of  incidence  up  to  30*  from  broadside  and  bistatio 
angles  up  to  15*.  Finally,  the  accuracy  of  Kell's 
relationship  proved  to  be  polarization  dependent  when 
surface  discontinuities  in  the  form  of  90  \  wedges  were 
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VERIFICATION'  AND  LIMITATIONS  OF  THE  MONOSTATTC-BTSTATIC 
RADAR  CROSS  SECTION  RELATIONSHIP  DERIVED  BY  KELL 


T  .  T ntrnduct i on 


Overvi pw 

One  method  used  t.o  reduce  the  detectability  of  an 
aircraft  by  radar  is  to  shape  it  ao  the  energy  from  the 
radar  is  reflected  away  from  the  radar’s  receiver.  By 
reflecting  the  incident  energy  away  from  the  receiver  the 
aircraft  presents  a  smaller  radar  cross  section  (RCS)  and, 
therefore,  is  more  difficult  to  detect.  Monostatic  radar,  a 
radar  in  which  the  radar’s  transmitter  and  receiver  antennas 
are  collocated,  is  more  common  than  bi static  radar,  a  radar 
with  separated  transmitter  and  receiver  antennas.  For  this 
reason,  most  aircraft  designed  to  present  a  low  RCS  reflect 
energy  away  from  the  direction  of  the  transmitting  source. 
However,  reducing  the  RCS  in  one  direction  using  shaping 
techniques  results  in  an  increased  RCS  in  another  direction 
(Knott,  1985:10).  Thus,  even  today’s  most  stealthy 
aircraft,  those  designed  to  hinder  detection  by  radar  at 
common  flight  aspects,  may  still  be  detectable  by  bistatic 
radar  receivers  (Adams,  1988:27). 

For  military  planners  to  project  a  weapon  system’s 
ability  to  penetrate  enemy  air  defenses,  it  is  necessary  to 
perform  both  monostatic  and  bi static  RCS  analyses. 
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Unfortunately,  however,  t.he  evolution  of  bistat.io  RCS 
prediction  and  measurement,  techniques  has  not  kept  pace  with 
that,  of  monostatic  techniques  (Hunka,  1978:243),  Because  of. 
the  threat  hi  static  radars  pose  for  even  stealth  aircraft 
and  the  need  to  plan  for  and  reduce  that  threat,  methods 
need  to  be  developed  that,  accurately  predict  bistat.io  RCS. 

Background 

The  earliest  radars  were  hi  static,  but  with  the 
appearance  of  monost.atic  radar  in  the  1930’s  bistatic  radar 
faded  from  prominence.  Bi static  radar  interest  surfaced 
again  in  the  mid  1950’s  leading  to  bistatic  radar 
applications  in  satellite  tracking  systems  and  guidance 
control  in  low  altitude  cruise  missiles  (Skolnik,  1961:20; 
Biggs  and  McMillen,  1979:1).  However,  the  development  of 
bist.atic  RCS  prediction  and  measurement  techniques  did  not 
receive  the  attention  monostatic  RCS  techniques  enjoyed 
until  later. 

Bistat.io  RCS  measurements.  Bi  static  RCS  measurements 
are  inherently  more  difficult  and  complex  than  monostatic 
measurements  simply  because  the  receiving  antenna  must  be 
moved  for  each  measurement  ( Hunka , 1 977 : 243 ) .  The  Air  Force 
Radar  Target  Scatter  Facility  (RATSCAT),  an  outdoor  radar 
range  at  Holloman  AFB,  New  Mexico,  uses  a  bistatic  mobile 
receiver  van  to  enhance  the  mobility  of  their  receiving 
antenna  (Dynalectron  Corp. ,  1985:43).  However,  they  must 
still  deal  with  the  complexities  involved  in  bistatic;  RCS 
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mensurement.fi  not.  suffered  by  monostatic  measurements.  Each 
repositioning  of  the  receiver  calls  for  recal ibrnt i on .  This 
means  the  target  of  interest  must,  be  replaced  with  a 
calibration  target  so  that  background  returns  can  be 
identified.  Multipath  scattering  must  also  be  accounted  for 
with  each  position  change.  And,  when  all  is  ready,  the 
target,  sometimes  as  large  as  a  full  scale  model  of  a  R1-R, 
must  be  remounted.  All  of  this  takes  a  great  deal  more  time 
and  costs  considerably  more  than  monost.atic  measurements 
( Tomko ,  1 988  )  , 

Ristatic  RCS  nredintion.  Before  the  1950’s,  the 
scientific  community  in  general  felt  the  theoretical 
calculation  of  a  complex  target’s  RCS  was  beyond  the 
capabilities  of  the  known  technology.  By  the  early  1950's, 
approximation  techniques  began  evolving  that  allowed  the 
estimation  of  the  RCS  for  such  complex  shapes  as  airplanes, 
missiles,  and  satellites.  Through  experience  and  the 
appearance  of  high  speed  computers  the  theoretical 
estimation  of  an  average  monost.atic  cross  section  circa  1958 
was  within  4  dB  of  measured  values  (Crispin  and  others, 
1968:v).  Bistatic  prediction  methods  were  also  developed 
during  this  period,  but  there  was  no  experimental  data  with 
which  to  compare  their  accuracy  (Siegel  and  others, 

1955:305) . 

One  could  easily  speculate  that  the  previously 


mentioned  problems  regarding  bi static  RCS  measurements  were 
the  reason  why  bi  static  RCS  experimental  data  was  in  short, 
supply.  Because  monostatin  RCS  measurements  are  much  easier 
to  perform  than  bi static  measurements,  if  a  means  of 
predicting  bi  static  RCS  from  monost.atic  measurements 
were  possible,  it  would  be  less  costly  and  faster  than 
performing  bi static  RCS  measurements.  This  approach  to 
bi static  RCS  prediction  would  be  beneficial  to  operations 
such  as  the  RATSCAT. 

Statement  of  the  Problem 

Tn  1965,  Robert  Kell  presented  a  relationship  between 
monost.atic  and  bistatic  radar  cross  section  which  appears  to 
offer  a  means  of  predicting  bistatic  RCS  using  monostatic 
measurements  under  certain  conditions.  The  purpose  of  this 
thesis  is  to  verify  the  theoretical  accuracy  of  Kell’s 
relationship  and  to  determine  the  physical  limitations  with 
which  Kell’s  method  has  practical  application. 

Scope 

The  basis  for  Kell’s  method  is  developed  first.  This 
development  begins  with  a  description  of  the  monost.atic- 
bist.atic  equivalence  theorem  (MBET),  the  foundation  of 
Kell’s  method.  Following  this  is  the  mathematical 
formulation  of  Kell’s  hypothesis. 

The  results  from  an  investigation  in  which  Kell’e 
method  is  theoretically  applied  to  specific  shapes  whose 
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bistat.ic  cross  section  can  be  accurately  modeled  are 
presented.  The  accuracy  of  Kell’s  method  when  subjected  to 
changes  in  either  bist.atio  angle,  electrical  size,  or 
surface  discontinuity  will  then  be  discussed. 

Approach 

The  limitations  of  Kell’s  method  can  be  f ound 
theoretically  by  comparing  computer  generated  bistatic  radar 
cross  sections  to  the  equivalent,  monost.atic  cross  sections 
produced  using  Kell’s  hypothesis.  Computer  programs  either 
in  use  by  Air  Force  agencies  or  developed  from  proven 
monost.atic  and  bi  static  RCS  models  will  be  used. 

Kell  suggests  his  method  is  most  accurate  in  the  high 
frequency  regime  t.o  about  10*  of  bi  static  angle  using 
targets  dominated  by  specular  returns  (Kell,  1965:987). 

This  investigation  will  apply  Kell’s  method  to  targets 
meeting  these  conditions  in  order  to  verify  the  accuracy  of 
this  method.  The  parameters  of  bistatic  angle,  electrical 
size,  angl ?  of  incidence,  and  surface  continuity  will  then 
be  varied  so  some  of  the  limitations  of  Kell’s  method  might 
be  brought,  out.  In  the  case  of  surface  continuity,  a  target 
with  90*  wedges  will  be  used  to  determine  the  limitations  of 
Kell’s  relationship  when  surface  discontinuities  are 
present . 
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Summary  of  Remaining  Chanters 


Chapter  II  is  a  development  of  Kell's  method.  Tt 
begins  with  a  description  of  the  monostatio-bistatie 
equivalence  theorem  (MRET).  A  review  of  the  literature 
regarding  the  MRET  is  then  provided.  This  is  followed  by 
the  mathematical  formulation  of  Kell's  method. 

Chapter  ITI  describes  RCS  models  for  the  sphere,  flat 
plate,  and  cylinder  for  both  the  bistatic  case  and  Kell's 
equi  valent  monost.atic  case.  Chapter  ITI  then  presents  the 
results  from  the  comparison  of  the  two  cases. 

Chapter  TV  contains  the  conclusions  reached  from  this 
study  and  recommendations  for  further  research  in  this  area. 
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T  T .  Theory 


Background 

Robert  Kell  implies  his  method  extends  the  monostatic 
bistat.io  equivalence  theorem  (MRET)  to  include  more  general 
oases  (Kell,  1965:983).  In  an  effort  to  offer  a  more 
complete  development  of  Kell’s  method,  this  chapter  first 
presents  the  MBET  and  documentation  supporting  the  concept 
of  predicting  hi  static  RCS  from  monostatic  RCS . 
Unfortunately,  there  has  been  little  research  performed  in 
this  area.  Tt  will  be  seen  that  the  MBET  relies  \ipon 
physical  optics  approximations,  and  therefore,  it  is 
applicable  to  only  those  targets  fitting  into  the 
limitations  set  by  physical  optics  techniques.  The 
development  of  Kell’s  method  shows  his  approach  relies  upon 
the  interaction  between  individual  scattering  centers.  Kell 
■ises  the  concept  of  reradiation  lobe  patterns  of  the 
individual  scattering  centers  to  define  the  bistatic 
scattering  pattern  in  terms  of  the  monostatic  pattern  and 
bi static  angle  (Kell,  1965:983). 

Monostatic-Bistatic  Equivalence  Theorem 

In  general,  discussions  of  bistatic  RCS  are  divided 
into  t.wo  areas.  The  first  is  concerned  with  bistatic  angles 
less  than  180*.  The  second  case  takes  intc  account  bist.atic 
angles  approximately  equal  to  180*.  This  is  known  as 


7 


forward  scattering.  Tt  has  significantly  different  behavior 
than  biatatio  scattering  and  is  not  addressed  in  this  study. 

In  considering  the  case  where  the  bistatic  angle  is 
less  than  1B0*,  Crispin  and  others  developed  what  has  come 
to  be  known  as  the  monos tat i e-bi static  equivalence  theorem 
(MRFT).  Referring  to  Figure  1,  this  theorem  claims  as  the 
wavelength  becomes  small  compared  to  target  sise  the 
bist.atio  cross  section  obtained  with  a  transmitter  direction  R 
and  a  receiver  direction  approaches  the  RCS  obtained  when 
the  transmitter  and  receiver  are  collocated  at  R  ♦ 
assuming  the  target  is  sufficiently  smooth  (Crispin  and 
others,  196R:1S8).  More  succinctly,  the  MBET  is  saying  the 
bi static  RCS  is  equal  to  the  monostatic  RCS  in  the  direction 
of  the  bisector  of  the  bistatic  angle  as  the  wavelength 
becomes  very  small. 


Bistatic-Monostatic  Relationship 
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Proof .  Crispin  offers  the  following  proof  for  the  MftRT 


(Crispin  and  others,  lPfift:lS8).  Radar  Croa*  Section  as  a 
function  of  transmitter  and  receiver  position  is  derived  as 
follows  (Siegel  and  others,  1955:10).  Assuming  the  surface 
of  the  snatt.ering  body  is  perfectly  conducting  the  equation 
for  the  scattered  magnetic  field  is  given  by 


H  - 

“  4  n 


x  Ht) 


x 


-ikR 

7  (  — - )ds 


(1  ) 


where  HRO  =  scattered  magnetic  field 

fi  =  the  unit  vector  normal  to  the  surface 

Ht  =  tangential  component  of  the  magnetic  field 
on  the  target’s  surface 

R  =  distance  between  the  receiver  and  the 
integral  on  point 

k  =  2n/x  (wave  number) 

S  ~  region  of  integration  (target’s  surface) 

By  assuming  the  incident  field  is  an  infinite  plane  wave, 
can  be  approximated  as  twice  the  tangential  component  of  the 
incident  magnetic  field  on  the  illuminated  side.  Letting 
the  incident  magnetic  field  have  a  magnitude  HQ  and 
direction  £  allows  Ht  to  be  written  as 

=  2itHne~^^  r  )  on  the  target’s  (2) 

illuminated  side 


=  0 


on  the  target’s  shadow  side 
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where  R  s  a  unit  vector  directed  from  the  transmitter, 

assumed  to  be  in  the  fer  field,  to  the  origin  of 
the  coordinate  aye tern 

r '  a  radius  vector  from  origin  to  any  point  on  the 
target’s  surface 

tt  a  fi  -  (ft  .  *)A 


Using  far  field  Approximations 


7(?ikR)  a  -ik>?ikR  )fiQ 
R  R^“ 


(3a) 


and 


R  a  R'  -  Cosa( 


(3b) 


where  R'  a  distance  from  the  origin  to  the  receiver. 

fiQ  =  a  unit  vector  directed  from  the  receiver  to  the 
origin 

fiQ  *  r' 

Coso  s  -  — — — 

lF*l 

Letting  HQ  =  1 ,  substituting  eqs.(2)  and  (3)  into  (1)  gives 


-ikR ' 


^sn  “  r-*  ^ 

(4) 

where 

F  =  ||[<n0  •  &)f  -  •  f)&] 

(5a) 

f  ikr  •  (fi0+R ) 

and  f  =1  fie  ds  (5b) 


a'  =  illuminated  region  of  the  target 
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Thr  RTS  can  now  be  expressed  as  (Crispin  and  others* 
1068:138) 

o  *  4*t  |FX|2  ♦  |FV|2  ♦  |Fj2  1  (6) 

In  Figure  1,  the  transmitter  is  placed  on  the  a  axis 
and  the  receiver  is  constrained  to  He  in  the  y-r  plane. 
With  this  configuration  the  monnat.atin  case  is  defined  by 

fi0  s  fSinS  -  SCosS  (?a) 


R  s  £ Sin#  -  ftCos#  (?b) 

A  s  RCose  ♦  ^SinaCosm  -t  RSinaSin*  (7c) 

whi  le  in  the  hist  at ic  case 


ftn  s  £Sin2*  -  6Cos2»  (8a) 

R  =  -R  (8b) 

&  =  RCosa  ♦  £Sina  (8c) 

Considering  the  vector  f  as  the  wavelength  becomes  small 
compared  to  the  target 


f 


I  fi0+R  I  V 


(9) 


For  the  monostatic  case,  eq.  (9)  becomes 


f 


(10) 
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And  in  the  bi static  rasp,  eq.  (9)  becoaes 

t .  *  f  d-ii,f'  •  ,*»K:n,*irt.  nn 

1  'a* 

where  fk  =  $Sine  -  ftCoae 

Evaluating  eqs. (10)  end  (11)  by  stationary  phene  yields 

f  *  {£ASine  *  f.ACoae]  exp(ikC)  (aonoatatic  case)  (12) 

f  2  t^ATane  -  ftA]  exp(ikr  Cnee)  (bistatic  case)  (13) 

Substituting  eqs.  (?)  end  (10)  for  the  aonoatatic  ease  into 
eq,  (5) 

F  =  YZ  [Aeikc($Cosa  ♦  ^SiaeCose  *  ISineSine]  (14) 

Applying  eq.  (14)  to  eq.  (6)  gives  the  aonoatatic  RCS 

om  =(kA)2/*  (15) 

Now,  substituting  eqs. (8)  end  (11)  for  the  bistetic  case 
into  eq.  (5) 

F  =  ^  (A  exp(ikC  Cose) ] [ ( SineSin2e) ( $Tane  -  2) 

-  (TaneSin2e  +  Cos2e)(ftCose  ♦  ^Sine)]  (16) 

Again,  applying  eq.  (16)  to  (6)  gives  the  bistatic  RCS 

ob  =  (kA)2/*  (17) 
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Comparing  eqs.  (15)  and  (17),  il  can  be  seen  that  aa 
wavelength  becomes  small  compared  to  target  sise  the 
monostatic  RCS  taken  on  the  bisector  of  the  hiatatic  angle 
is  equivalent  to  the  hiatatic  RCS. 

The  equations  the  MRET  is  developed  from,  eqs.  (1)  and 
(2),  are  derived  using  physical  optica  techniques  (Siegel 
and  others,  1953:10).  Therefore,  the  MBFT  is  valid  only  if 
the  restrictions  required  by  physical  optics  are  satisfied. 
Those  restrictions  are  given  below  (Skinner  and  Jost, 

1 988 : 2 ) : 

1.  The  target  must  have  a  large  enough  RCS  that  higher 
order  effects  do  not  make  a  significant  difference. 

2.  The  target  must  be  significantly  larger  than  a 
wavel  engt.h . 

3.  A  significant,  amount  of  the  target's  surface  should 
be  within  approximately  45*  of  the  normal  of  the  bisector  of 
the  transmitter  and  receiver. 

4.  The  target  must  be  perfectly  conducting. 

These  restriction  limit  the  applicability  of  the  MBET.  Two 
studies  were  found,  however,  in  which  the  MBET  was  used  to 
determine  bistat.ic  RCS. 

Documentation  supporting  the  MBET.  Biggs  and  McMillen, 
in  their  efforts  to  formulate  the  bistatic  RCS  of  a  prolate 
spheroid,  found  the  MBFT  returned  accurate  results.  The  two 
discrepancies  these  researchers  noted  when  comparing 
measured  data  with  data  obtained  using  the  MBET  were  within 
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one  percent,  and  were  attributed  to  the  fact  the  MBFT 
requires  the  target  be  infinitely  removed  from  the 
transmitter  and  the  receiver.  When  it  is  not,  the  hi  static 
receiver  is  not,  looking  at  the  exact  same  spot  as  the 
monost.atio  receiver  (Riggs  and  McMillen,  1979:1  3). 

Coleman,  while  testing  an  aircraft.  RCS  computer 
prediction  method,  MTSCAT  ITT, -against  measured  data  of  the 
F-5  aircraft,  found  hi  static  angles  up  to  16*  had  little 
effect  on  the  RCS.  This  made  the  MRET  attractive.  By 
taking  advantage  of  this  approximation  he  was  able  to  reduce 
program  running  time  for  the  bi static  case  for  small 
bist.atio  angle  to  that  of  the  inonostatic  case.  At  the  same 
time,  he  was  able  to  resolve  some  inaccuracies  by  avoiding 
certain  elements  within  the  bistatic  RCS  derivation 
(Coleman,  1977:34). 

Tn  the  process  of  applying  the  MBET,  Coleman  developed 
error  estimates  for  parallel,  perpendicular,  and  cross 
polarizations.  His  interpretation  of  the  error  estimates 
reach  the  same  conclusions  that.  Kell  does  in  the  development 
of  his  approach.  First,  the  RCS  is  determined  by  the 
relationships  between  scattering  center.  The  second,  that 
the  equivalent  monostatic  RCS  includes  the  factor  Cos((J/2) 
where  8  is  the  bistatic  angle  (Coleman,  1977:43). 

Kell's  Method 

Kell  developed  a  method  to  obtain  bistatic  RCS  using 
monostatic  RCS  measurements.  His  method  uses  the  concept  of 
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reradiation  lobe  patterns  from  individual  scattering  centers 
to  define  the  bi static  RCS  in  terms  of  monostatic  patterns 
and  the  bi static  angle.  Kell’s  approach  uses  the  MBET  by 
relating  the  bi  static  RCS  to  the  monostatic  RCS  viewed  on 
the  bisector  of  the  bistatie  angle  by  a  factor  of  Cos(0/2) 
(Kell,  1965:983). 

Mathematical  formulation.  The  geometry  Kell  used  is 
given  in  Figure  2. 


Figure  2.  Bistatic  Coordinate  System  with  Antennas 
in  the  x-z  Plane 

It  oan  be  seen  in  Figure  2  that  Kell  aligns  the  bisector  of 
the  the  bistatie  angle,  t,  and  the  centroid  of  the  target 
along  the  z  axis.  Applying  the  Stratton-Chu  equation  for 
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the  scattered  magnetic  field  to  this  geometry,  Kell  arrives 
at  the  following  expression  for  the  bistatic  RCS  (Kell, 
1965:985) 


o 


n 


exp  [i  k2zCos  (  0/2  )]  dz 


(  18a) 


where 


T  (z)  =  |[(n  x  hs)  x  fD(e  0  )  +  (fi»hs)  (ei0  )  f  „ 

-  (n  x  egje1^]  x  p(©,z)  de 


(18b; 


The  derivation  of  eq.  (18)  can  be  found  in  appendix  A. 

Kell  made  the  following  observations  regarding  eq.  (18) 
(Kell,  1965:985): 

a)  I(z)  is  an  exact  description  of  the  surface 
geometry,  illuminating  and  observing  ray  geometry,  and 
surface  wave  propagation  effects 

b)  While  eq.  (18)  has  the  form  of  a  physical  optics 
expression  it  is  not  a  physical  optics  expression  because 
T(z),  in  this  case,  is  exact  as  opposed  to  a  physical  optics 
approximation 

c)  Eq.  (18)  can  be  divided  into  a  sum  of  subintegrals 
whose  end  points  correspond  to  scattering  centers.  The  end 
points  of  each  subintegral  are  determined  by  the  range  in  z 
over  which  the  respective  integrands  are  continuous. 
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According  to  Kell,  the  phaRe  of  a  given  scattering 
center  plays  a  large  role  in  the  contribution  that 
scattering  center  makes  to  the  total  RCS .  The  magnitude  of 
the  contri butr i on  of  an  individual  scattering  center  to  the 
overall  return  is  directly  related  to  the  area  covered  by 
the  scattering  center.  The  size  of  the  area  of  the 
scattering  center  is  directly  related,  in  turn,  to  the  phase 
of  the  scattering  center.  All  points  with  phases  within  n/2 
of  the  scattering  center’s  phase  contribute  positively  with 
the  scattering  center’s  return.  The  greater  the  area  within 
a/2  of  the  scattering  center’s  phase  the  greater  the 
contribution  from  that  scattering  center.  This  concept 
applies  equally  t,o  both  monostatic  and  bi static  RCS 
(Kell,  1965:985). 

Differences  will  often  arise  between  bistat.ic  RCS  and 
mnnostat.io  RCS  viewed  on  the  bisector  of  the  bistatic  angle 
as  the  bistatic  angle  changes.  Kell  relates  these 
differences  to  the  following  changes  in  the  scattering 
centers  (Kell,  1965:985): 

a)  Changes  in  the  phase  of  scattering  centers  relative 
to  other  centers 

b)  Changes  in  strength  of  the  radiation  from  a  given 
scattering  center 

c)  Change  in  the  number  of  scattering  centers. 

Looking  at  eq.  (18),  phase  is  clearly  dependent  on  the 

biBtatic  angle  by  a  factor  of  Cos(0/2).  Kell  couples  this 
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factor  with  the  wave  number  k  and  then  views  a  change  in 
bi static  angle  as  a  change  in  wavelength.  As  the  relative 
wavelength  changes  with  bistatic  angle  the  area  within  n/2 
wavelengths  of  the  scattering  centers  also  changes  (Kell, 
1965:985) . 


The  radiation  strength  of  a  scattering  center  is 
dependent  on  the  strength  of  the  surface  current,  among 
other  things,  at  that  location.  With  changes  in  the 
bi static  angle  there  can  be  changes  in  the  tangential 
components  of  the  incident  field  resulting  in  a  change  n 
the  surface  currents. 

Kell  describes  the  appearance  and  disappearance  of 
scattering  centers  as  aspect  angle  and  target  geometry 
dependent,  and  declines  further  explanation  (Kell,  1965:985). 

The  determination  of  bistatic  RCS .  Kell  approximates 
the  target  by  breaking  it  into  a  collection  of  discrete 
scattering  centers.  Those  centers  directly  illuminated  by 
the  incident  wave  are  termed  simple  centers  while  those 
illuminated  by  reflections  from  other  parts  of  the  target 
are  termed  reflex  centers.  Kell’s  addresses  only  the  simple 
scattering  centers.  Using  the  concept  of  discrete 
scattering  centers,  Kell  then  derives  the  total  RCS  by 
summing  the  product  of  the  individual  center’s  RCS  and  a 
relative  phase  term  (Kell,  1965:986): 


a  = 


M 

L  !<o m) 

m=  1 


1/2 


exp(  i  ^jp)  ] 


(19) 
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where  o  =  total  RCS 

M  =  number  of  discrete  scat.terers 

om  =  RCS  of  the  m*'*1  scattering  center 

=  phase  of  field  scattered  by  the  center 

relative  to  that  scattered  by  the  first  center, 
and  where  om  is  evaluated  for  the  desired 
polari zation . 

Tn  examining  the  relations  between  the  monost.atic  and 
hi  static  RCS  for  a  given  aspect  angle  ©,  Kell  expresses  the 
individual  scattering  center’s  phase  in  eq.  (19)  as  a  t.wo- 
term  sum  (Kell,  1965:986): 

om  =  2kzm( ®)Cos ( 8/2 )  +  Km  (20) 

where  zm(e)  =  distance  between  the  m^  and  the  first  phase 

center,  projected  on  the  bistatic  bisector 
ax*  is 

=  residual  phase  contributions  of  the  m**^ 
center 

The  hi  static  RCS  for  scattering  centers  that  do  not  move 
with  changes  in  bistatic  angle  can  now  be  written  as  (Kell, 
1965:986) 

M  2 

0(0)  =  £  (om)1//2  exp[  i2kzn,(  ©)Cos  (  fi/2 )  +  iCm]  (21) 

m=  1 

Kell  uses  eq,  (21)  to  make  a  statement  for  monostatic- 
bistatic  equivalence  and  the  conditions  required  for  its 
occurrence. 

If  for  a  chosen  aspect  angle,  the  following  conditions 

hold : 
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1)  The  RCS  may  be  written  as  a  squared  sum  of  fields 
from  discrete  scattering  centers,  and 

2)  The  amplitude  om  ,  position  zm,  and  residual  phase 
Cm  are  insensitive  to  the  bist.atir  angle  0  over  the 
range  of  0  considered,  for  those  centers  which  are 
significant,  members  in  this  sum; 

it  then  follows  that  the  bi static  cross  section  of 
aspect  angle  ©  and  bistat.ie  angle  B  is  equal  t.o  the 
monostatic  cross  section  measured  on  the  bisector  at  a 
frequency  by  the  factor  Cos (0/2)  [Kell,  1965:98?]. 


Implementing  Kell’s  method.  The  st.ep-by-step  process 
to  employ  Kell’s  method  calls  for  first  measuring  the 
monostat.io  cross  section  as  a  function  of  aspect  angle  and 
at  a  frequency  higher  than  the  one  desired  for  the  bistat.ie 
data  by  a  factor  of  Secant.  ( 0/2 )  .  The  reason  Kell  raises  the 
frequency  can  be  seen  in  eq.  (18).  As  previously  explained, 
phase  changes  are  indicated  by  the  factor  Cos (0/2) 
accompanying  the  wave  number  k.  This  can  be  considered  a 
change  in  wavelength  or  frequency  corresponding  to  a  change 
in  bistat.ie  angle.  The  measured  monostatic  data  is  then 
translated  along  the  aspect  angle  axis  0/2  degrees  and  the 
measurement,  frequency  is  reduced  by  Cos  (0/2)  (Kell, 

1965:987)  . 


Summary 

It  was  demonstrated  theoretically  that  the  MBET 
provides  a  relationship  between  monostatic  and  bistatic  RCS 
when  physical  optics  approximations  are  valid.  The  utility 
of  the  MBET  for  small  bistatic  angles  was  shown  in  the  work 
performed  by  Coleman. 
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Kell’s  method  approaches  the  problem  of  equaling 
bistatic  RCS  and  monostatic  RCS  in  much  the  same  way  as  the 
MBFT  in  that  it  uses  the  same  scattering  geometry  and  is  in 
the  high  frequency  regime.  Unlike  the  MBET  which 
approximates  surface  currents  to  determine  RCS,  the 
theoretical  basis  for  Kell’s  method  lies  in  the  summation  of 
the  contributions  from  individual  scattering  centers.  His 
mathematical  formulation  of  this  hypothesis  is  exact  up  to 
the  point  where  he  makes  far  field  approximations  and 
evaluates  the  bist.atir  RCS  integral  using  stationary  phase 
mehtods.  According  to  Kell,  the  change  in  the  contributions 
from  the  scattering  centers  with  changing  bistatic  angle  is 
due  in  large  part  to  changes  in  the  relative  phase  of  each 
scattering  center.  For  small  bistatic  angles,  Kell  reports 
phase  will  change  by  a  factor  of  Cos (0/2).  Kell  noticed 
that  the  Cos(/9/2)  factor  occurs  with  the  wave  number  k. 
Viewing  this  as  a  change  in  freqxiency,  Kell  hypothesizes 
that  monos tat ic  RCS  is  equal  to  bi static  RCS  when  the 
frequency  of  the  monostatic  case  is  increased  by  a  factor  of 
Sec ( 0/2 ) . 

Kell’s  method  could  be  tested  theoretically  by 
comparing  computer  predicted  monostatic  cross  sections  to 
bistatic  RCS  predictions. 


ITT.  Results 


There  are  proven  RCS  prediction  methods  that  can 
accurately  approximate  bistatic  RCS  as  well  as  monostatic 
RCS  for  certain  shapes.  Kell’s  hypothesis  can  be  tested  by 
comparing  data  derived  from  bist.at.ic  prediction  methods  with 
that  obtained  from  monostatic  approximations  at 
appropriately  adjusted  frequencies  and  aspect  angles.  Three 
simple  shapes  are  used  to  prove  the  accuracy  and  limitations 
of  Kell’s  hypothesis.  These  are  the  sphere,  square  flat 
plate,  and  right  circular  cylinder.  An  exact  solution  is 
used  to  predict  the  RCS  of  a  sphere,  physical  optics  methods 
are  applied  to  approximate  the  RCS  of  a  square  flat  plate, 
and,  finally,  the  RCS  of  a  right  circular  cylinder  is 
obtained  using  geometrical  theory  of  diffraction  techniques. 

Testing  Kell’s  Hypothesis 

Kell’s  hypothesis  can  be  tested  theoretically  by 
comparing  the  predicted  bistat.ic  RCS  to  the  predicted 
monostatic  cross  section  at  the  bisector  of  the  bistatic 
angle.  The  bistatic  RCS  is  calculated  for  each  target  at 
designated  angles  of  incidence  with  the  bistatic  angle 
varying  from  0*  to  60*  in  1*  increments  for  the  square  flat 
plate  and  the  cylinder  while  the  range  of  bistatic  angles 
for  the  sphere  is  0*  to  90*  in  1*  increments.  All  bistatic 
RCS  predictions  are  made  at  a  constant  frequency. 

Monostatic  RCS  predictions  are  made  at  aspect  angles  equal 
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to  the  designated  angle  of  incidence  plus  half  the  bist.atie 
angle.  The  monost.atic  RCS  prediction  frequency  is  increased 
by  a  factor  of  secant (0/2)  (B  =  bi static  angle)  above  the 
frequency  used  for  the  bi static  case. 

The  accuracy  of  Kell's  theory  can  be  determined  by 
comparing  the  data  obtained  from  each  method.  A  plot  of 
bi  static  RCS  and  monostat.ic  RCS  versus  bistatic  angle  at  a 
given  angle  of  incidence  shows  how  closely  Kell's  equivalent 
monost.atio  RCS  coincides  with  the  bistatic  RCS  as  bist.atic 
angle  increases. 


Sphere 

Approach .  In  the  case  of  a  perfectly  conducting 
sphere,  the  exact  solution  for  the  bi static  cross  section  of 
a  plane  electromagnetic  wave  is  taken  from  the  Mie  series. 

It  is  (Bowman,  1969:400) 


o(e,«)  =  [  |Sj  (e) !  2Cos2(  a)  +  |s2(e)|2Sin2U)]  (22) 


where , 


S.  = 


.  v-  (-1)  2n+l  (\_  n' 
’  * 1  n  =  l  n^fT+il  I  n  ae 


ap1 (Cose) 


p1 (Cose) 

an  Sin.  ]  (23> 


„  _  .  V  ( -1 )n ( 2n  +  l )  f 

2  "  nTl - n(nTiribn^I 


P1 (Cose) 


ne 


—  a 


3P1 (Cose) 
n 


n  f  e 


]  (24) 
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•  *  bistatic  angle 


+  *  angle  between  plane  of  scattering  and  the  plane 
containing  the  incident  i  field  and  the  direction 
of  incidence 

and  (Knott,  1985:90) 


jn(k») 


(25) 


b  = 
n 


kaj  . (ka)  -  nj  (ka) 

H"  1  n 

kah (  ( ka )  -  nhf  ^ ( ka ) 

n  - 1  n 


(26) 


wi  th 


jn(ka)  =  spherical  Bessel  function 

h^^(ka)=  spherical  Hankel  function  of  the  first  kind 

k  =  2*/x  (the  wave  number) 
a  =  sphere  radius 

For  the  monostatic  case,  eq,  (22)  reduces  to  (Bowman, 
1969:401 ) 


_  _  4a  ^  ( 1-  )n2n+l  _  ,1 

'  n=  1  n(B+1)  n> 


(27) 


The  geometry  describing  scattering  from  a  sphere  is 
given  in  Figure  3. 
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Figure  3.  Scattering  Geometry  for  a  Sphere 

(Knott*  1985:89) 


Two  computer  programs  given  in  Appendix  B  were  used  to 
derive  the  RCS  of  various  sixed  spheres.  One  progrsm*  a 
product  of  Ohio  State  University  that  is  used  on  the  Wright- 
Patterson  AFB  indoor  radar  range  for  calibration  purposes, 
gives  an  exact  solution  for  the  radar  cross  section  of  a 
sphere  for  any  number  of 

bistatxc  angles.  For  this  study,  the  bistatic  angle  was 
varied  from  0*  to  90*  in  increments  of  1*.  The  second 
program  is  a  modified  version  of  the  first  progrsm.  It 
computes  the  monostatic  RCS  of  s  sphere  over  s  range  of 
discrete  frequencies  with  each  frequency  increased  by  a 
factor  of  aecant(i/2)  over  the  previous  frequency.  8  is  the 
bistatic  angle  and  ranges  from  0*  to  90*  in  increments  of 


1*.  Varying  the  frequency  in  this  way  directly  relates  the 
equivalent  monnstatie  RCS  data  in  a  one-to-one 
correspondence  to  the  bi static  RCS  obtained  for  each 
biatstic  angle. 

Figures  4  through  27  show  the  relationship  between 
bi static  RCS  and  the  RCS  obtained  using  Kell’s  method.  Six 
sphere  sixes  of  ka  equal  to  1 ,  3,  8,  10,  30,  and  50  at  10 
GH*  were  exaained  using  horixontal  and  vertical 
polari eations . 

Near  resonant  sixe  sphere  results.  A  comparison  of 
Figures  5  and  ?  show  there  is  a  significant  difference 
between  vertically  and  horixontally  polarised  bistatic  RCS 
and  between  biatstic  and  aonostatic  RCS  for  an  electrically 
small  sphere  (ka  si).  The  differences  are  caused  by 
creeping  waves.  For  a  vertically  polarised  field  as  shown 
in  Figure  4,  boundary  conditions  for  a  perfectly  conducting 
surface  require  the  tangential  component  of  the  surface 
electric  field  to  be  sero  in  the  x-y  plane. 


Figure  4.  Creeping  Wave  Scattering  in  the  x-y  Plane 


26 


Therefore,  vertically  polarised  creeping  waves 
contribute  to  the  nonostatic  RC8,  but  aa  the  bist.atic  angle 
increases  in  the  x-y  plane,  creeping  waves  no  longer  scatter 
in  the  x-y  plane  and  the  bi static  RCS  decreases. 

For  horizontal  polarization,  Figure  7  deaonstrat.es  the 
creeping  wave  continues  to  contribute  as  the  bist.atic  angle 
increases.  In  the  case  of  ka  s  1  (Figure  7),  the  phase  of 
the  creeping  wave  is  relatively  insensitive  to  changes  in 
histat.ic  angle.  This  allows  the  creeping  wave  contribution 
to  equal  that  of  the  monostatic  case  at  bi static  angles 
greater  than  45*. 

The  effects  of  creeping  waves  can  also  be  seen  in  the 
RCS  data  for  spheres  ka  s  3  (Figures  9  and  11)  and  /ta  s  5 
(Figures  13  and  15).  The  horizontally  polarized  creeping 
wave  continues  to  contribute  to  the  RCS,  however,  oreeping 
wave  phase  is  becoming  more  sensitive  to  changes  in  bistatic 
angle  as  the  sphere  becomes  larger.  Bi static  and  monostatic 
RCS  are  in  agreement  to  approximately  35*  of  bistatic  angle 
for  ka  s  3  and  only  8*  of  bistatic  angle  for  ka  =  5.  For 
the  case  of  vertical  polarization,  the  bistatic  and 
monostatic  RCS  begin  to  diverge  at  1  to  2  degrees  of 
bi static  angle  for  ka=3,5. 

This  data  demonstrates  Kell's  relationship  cannot 
account  for  creeping  waves  that  are  polarization  dependent 


Failure  5.  Sphere  (ka  =  1;  Vertical  Polarization): 

Bi static  RCS  and  Equivalent  Monoatatic  RCS 


Figure  6.  Sphere  (ka=l;  Vertioal  Polarization): 

Difference  in  RCS  (Biatatic  -  Monoatatio) 
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Figure  7.  Sphere  (kasl;  Horieontal  Polarisation): 

Biatatic  RCS  and  Equivalent  Monostatic  RCS 


Figure  8.  Sphere  (kasl;  Horisontal  Polarisation): 

Difference  in  RCS  (Bistatic  -  Monostatio) 
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Figure  9.  Sphere  (ka=3;  Vertical  Polarisation#: 

Bistatic  RCS  and  Equivalent  Monostatic  RC8 


Figure  10.  Sphere  (ka=3;  Vertioal  Polarisation): 

Difference  in  RCS  (Bistatio  -  Monoatatic) 
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at  static 


at  static  amu  COTA) 


Figure  11.  Sphere  (ka:3:  Horizontal  Polarization): 

Bistatic  RCS  and  Equivalent  Monoatatio  RCS 
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Figure  12.  Sphere  (kas3;  Horizontal  Polarization): 

Pifferenoe  in  RCS  (Bintatio  -  Monostatic) 
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Figure  13.  Sphere  (kasS;  Vertical  Polarisation): 

Biatatic  RCS  and  Equivalent  Monostatic  RCS 


Figure  14.  Sphere  (ka&5;  Vertioal  Polarisation): 

Difference  in  RCS  (Biatatio  -  Monostatio) 
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Figure  15.  Sphere  (ka=5;  Horizontal  Polarization): 

Bistatic  RCS  and  Equivalent  Monostatic  RCS 
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Large  snhere  results.  The  equivalent  monostatic  cross 
section  values  are  very  close  to  those  predicted  for  the 
bistatic  case.  Much  of  the  slight  differences  between  the 
two  cases  are  caused  by  creeping  waves.  This  can  be  seen  in 
Figures  17,  19,  21,  23,  25,  and  27.  Attenuation  due  to  the 
larger  electrical  path  lengths  of  these  larger  spheres 
causes  the  reduction  in  creeping  wave  effects.  Looking  at. 
Figures  18,  20,  22,  24,  26,  and  28,  Kell’s  method  is  very 
good  well  beyond  10*  for  doubly  curved  surfaces  with  no 
surface  discontinuities  and  at  high  frequencies. 
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Figure  18«  8phere  (ka=10;  Vertical  Polarisation): 

Difference  in  RCS  (Bi static  -  Monostatio) 
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Figure  20.  8phere  (kaslO;  Horisontal  Polarisation): 

Difference  in  RCS  (Biatatic  -  Monostatic) 
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Figure  21.  Sphere  (kas30;  Vertical  Polarisation): 

Biotatic  RCS  and  Equivalent  Monoatatio  RCS 


Figure  22. 


8phere  (kas30;  Vertical  Polarisation): 
Difference  in  RCS  (Bistatio  -  Monoatatio) 
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Figure  23.  Sphere  (ka=30;  Horisontal  Polarisation): 

Bistatic  RCS  and  Equivalent  Monoatatio  RCS 


Figure  24.  8phere  (kasSO;  Horimontal  Polarisation): 

Diffaranoa  in  RCS  (Bistatio  -  Monoatatio) 
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Figure  25.  Sphere  (kasSO;  Vertical  Polarisation): 

Biatatio  RCS  and  Equivalent  Monostatic  RCS 
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Figure  28.  Sphere  (ka*50;  Horisontal  Polarisation): 

Difference  in  RCS  (Dietetic  -  Monoetetio) 
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E1&1  flax? 


Approach .  Using  the  geoaetry  given  in  Figure  29,  the 
hi  static  RCS  for  a  perfectly  conducting  square  flat  plate  la 
expressed  as 


_  _  *  _  ra  - -  Sint  (ha/2)  Sin* 

°  -  4"  l~  =»•<•♦»>  iiii/sHiir,*  ♦ 


♦  Sin(e+f ) ) 


2 

A] 


(28) 


where 


•  *  angle  of  incidence 
0  =  bistatic  angle 
a  =  length  of  aide 


Figure  29.  Biatatic  Scattering  Geoaetry  for  a 
Square  Flat  Plate 

The  equivalent  aonostatic  cross  section  uses  a 
different  geoaetryi  shown  in  Figure  30,  because  Kell* a 
aethod  requires  the  aonostatio  aspect  angle  bisect  the 
bistatic  angle.  This  changes  eq.  (28)  to 
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o  s 


(29) 


Figure  30.  Equivalent  Monostatic  Scattering  Geometry 
for  a  Square  Flat  Plate 

The  derivation  of  eqs.(28)  and  (29)  can  be  found  in 
appendix  C.  These  are  physical  optics  approximations. 

Using  physical  optics  to  compute  the  RCS  of  a  flat  plate 
limits  the  size  of  the  angle  of  incidence  that  can  be 
accurately  modeled  to  about  40*  from  broadside  (Knott, 
1985:174).  This  is  because  physical  optics  does  not  account 
for  edge  effects  (Skinner  and  Jost,  1988:8). 

Only  one  polarization  is  used,  E  field  in  the  x-z 
plane.  Beoause  edge  effects  are  not  considered  and  the 
plate  is  square  as  opposed  to  rectangular,  Kell’s  method  can 
be  adequately  described  with  one  polarization. 
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Flat  plate  results.  Tt.  can  be  seen  in  Figures  31 
through  42  that  Kell’s  method  degrades  as  the  angle  of 
incidence  moves  away  from  broadside.  Figure  31,  broadside 
incidence,  shows  Kell’s  method  holds  beyond  20*  of  bist.atio 
angle.  When  the  angle  of  incidence  is  between  8*  and  24* 
from  broadside,  Figures  33,  35,  and  37,  the  lobe  maximums 
are  in  close  agreement  to  about  16*  of  bistatio  angle.  The 
nulls,  however,  are  significantly  different.  The  large 
differences  in  nulls  cause  the  oscillation  between  0*  and 
20*  of  bistat.ic  angle  in  the  plots  showing  the  difference 
between  bistat.ic  and  equivalent  monost.at.ic  RCS,  Figures  34, 
36,  and  38. 

The  agreement  between  bi static  and  equivalent 
monostatic  RCS  breaks  down  before  10*  of  bistatic  angle  as 
the  angle  of  incidence  from  broadside  becomes  large.  From 
Figures  39  and  41,  significant  differences  begin  to  occur  at 
t.he  first,  null  as  was  seen  at  smaller  angles  of  incidence; 
the  second  lobes,  however,  do  not  agree. 
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Figure  31.  Square  Flat  Plate  (e  =  O'):  Bistatic  RCS 
and  Equivalent  Monostatic  RCS 


Figure  32*  Square  Flat  Plate  (•  s  0*) 

Difference  in  RCS  (Bistatic  -  Monoatatio) 
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Figure  33.  Square  Flat  plate  (•  =  8*):  Bistatic  RCS 
and  Equivalent  Monostatic  RCS 


Figure  34.  Square  Flat  Plate  (•  s  8*): 

Difference  in  RCS  (Biatatic 


-  Honoatatio) 
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Figure  35.  Square  Flat  Plate  (e  =  16*):  Bistatic  RCS 
and  Equivalent  Monostatic  RCS 
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Figure  36.  8quare  Flat  Plate  (•  s  16*): 

Difference  in  RCS  (Biatatic  -  Monoatatio) 
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Figure  37.  Square  Flat  Plate  (•  =  24*):  Bistatic  RCS 
and  Equivalent  Monostatic  RCS 


Fifura  40.  Squirt  Flat  Plata  (a  «  38*): 

Diffarenoa  In  RC8  (Blatatlo  -  Monoatatio) 
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Figure  41.  Square  Flat  Plate  (•  =  40*):  Bi static  RCS 
and  Equivalent  Monostatic  RCS 


Figure  42.  8quare  Flat  Plate  (•  s  40*): 

Difference  in  RC3  (Bi static  -  Monostatio) 
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Cvl i nrier 


Approach .  The  RCS  of  a  perfectly  conductinf  right 
circular  cylinder  ia  dominated  by  diffraction  at  moat  aapect 
anglea.  For  this  reason,  the  geometrical  theory  of 
diffraction  (GTD)  is  used  to  calculate  its  cross  section 
(Anderson,  1965:3-21).  According  to  GTD,  the  illuminated 
edges  of  the  cylinder  shown  in  Figure  43  as  Sj,  S2,  and  S3 
are  t.he  major  contributors  to  the  cylinder’s  RCS  (Anderson, 
1965:3-21). 


Figure  43.  Cylinder  Geometry 

The  directions  of  incidence  and  reflection  are  constrained 
to  lie  in  the  x-y  plane.  The  electrical  polarization  is 
vertical  when  in  the  z  direction  and  horizontal  when  in  the 
x-y  plane  (Anderson,  1965:3-21). 

The  RCS  of  a  perfectly  conducting  right  ciroular 
cylinder  is  expressed  in  eq.  (30)  (Anderson,  1965:3-22). 
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-  h(Cos«  +  Cos(e+jj) )  ] } 
x  { [Cos ( 2«/3 )  -  Cos(2fl/3) J'1 

T  [Cos  (  2  x/3  )  -  Cos(2(  n-2«-j5)/3)  ]-1} 

angle  of  incidence  in  degrees  from  end-on 

bistatic  angle 

half  the  cylinder’s  length 


a  =  cylinder’s  radius 

k  =  2«/x  (wave  number) 

Analytic  singularities  occur  at.  •  =  0*  and  0  s  90* 
in  eqs .  (31a)-(31c).  Therefore,  at  •  =  0*.  the  RCS  is 

computed  using  eq.  (32)  (Anderson,  1965:3-24). 

o  =  *[ka2Cos( 0/2) ]2  (32) 

At  ®  =  90*,  because  only  edges  Sj  and  S2  are  illuminated, 
eq.  (30)  reduces  to  (Anderson,  1965:3-24) 

o  =  4kah2Cos(fl/2)  (33) 

The  derivations  of  eqs.  (31a, c,b),  (32), and  (33)  are  given 
in  appendix  D.  Computer  data  generated  from  eqs.  (31 )— (33) 
Was  compared  to  measured  data  for  the  monostatic  and 
bi static  case  ( B  =  10*).  The  predicted  monstatic  case  was 
within  3-dB  of  the  measured  data  between  4*  and  86*  of  end- 
on  while  the  predicted  bistatic  case  was  within  5-dB  over 
the  same  range  of  aspect  angles. 

Cylinder  results.  It’s  obvious  something  is  amiss  when 
viewing  Figures  44  and  45  in  which  the  angle  of  incidence  is 
end-on.  The  disparity  arises  from  how  the  two  cross 
sections  are  calculated.  In  the  bistatic  case,  eq.  (32)  is 
used  for  all  bistatic  angles.  The  equivalent  monostatic 
case  calls  for  eq.  (32)  when  the  monostatic  angle  of 
incidence  is  equal  to  0*,  but  when  beta  does  not  equal  0*  it 
uses  eq.  30.  In  doing  this,  the  bistatic  RCS  is  calculated 
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as  if  the  target  is  s  disc  which  allows  the  use  of  the 
physical  optics  approximation  given  by  eq.  (32).  The 
equivalent  monostatic  RCS  is  calculated  the  same  way  when 
the  angle  of  incidence  is  end-on,  but  once  the  equivalent 
monostatic  angle  of  incidence  is  no  longer  exactly  end-on 
the  target  becomes  a  cylinder  and  GTD  approximations  are 
applied.  Tn  short,  scattering  for  the  equivalent  monostatic 
case  is  being  modeled  as  either  purely  specular  or  purely 
diffraction.  Neither  of  these  conditions  is  likely  at  small 
angles  from  end-on.  For  this  reason,  no  conclusions  can  be 
drawn  regarding  the  accuracy  of  Kell’s  method  when  either 
the  angle  of  incidence  or  bistatic  bisector  are  end-on. 
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Figure  44.  Cylinder  (•  *  0*; Vertical  Polarisation): 

Bi static  RCS  and  Equivalent  Monostatic  RCS 
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Figure  45. 


Cylinder  (•  *  0*»  Horisontal  Polarisation) 
Difference  in  RC8  (Bistatic  -  Monostatic) 
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It  can  be  aeen  in  Figures  47-65  that  Kell's  Method  ia 
polarisation  dependent.  Vertically  polarised  eroaa  sections 
are  eaooth  curves  while  horisontal  polarisation  produoea 
rapid  fluctuation  with  deep  nulla.  This  results  froa  the 
fart  that  the  scattered  fields  froa  a  right-circular 
cylinder  are  polarisation  dependent .This  is  deaonstrated  in 
Figure  46. 


uiiwtm  mrci  m*u  •a(ucnu) 

Figure  46.  GTD  Applied  to  Sight-Circular  Cylinder  to  Show 
Contributions  froa  Xlluainated  Sdges 

(Anderson*  1965:3-28) 
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The  plots  shown  in  Figure  46  depict  the  RCS  magnitude  of 
each  of  the  major  scattering  centers  with  respect  to  aspect 
angle.  The  vertically  polarized  RCS  is  dominated  by  the 
return  from  except  at  end-on  and  broadside.  Therefore, 
the  fluctuations  in  the  cylinder’s  vertically  polarized  RCS 
are  governed  by  changes  in  phase  from  Sj  as  the  aspect  angle 
changes.  This  results  in  relatively  smooth  curves  which  can 
be  seen  in  Figures  47,  51,  55,  59,  and  63. 

When  the  polarization  is  in  the  horizontal  direction, 
Figure  46  shows  that  the  RCS  is  no  longer  dominated  by  one 
scattering  center  but  is  a  combination  of  two  or  three 
scattering  centers  at  all  aspect  angles.  With  each  change 
in  aspect  angle  there  will  be  a  change  in  the  phase  of  the 
scattered  fields  from  each  scattering  center.  Combining  the 
: hases  from  each  scattering  center  produces  an  overall  phase 
that  is  sensitive  to  aspect  angle  changes.  This  results  in 
rapid  fluctuations  in  RCS  which  can  be  seen  in  Figures  49, 
53,  57,  61,  and  65.  Because  of  the  sensitivity  to  aspect 
angle,  the  equivalent  monostatic  case  does  not  track  the 
bistatic  case  well. 

Kell’s  relationship  holds  up  better  for  vertical 
polarization  because  the  relative  phase  between  the 
cylinder’s  scattering  centers  is  less  sensitive  to  changes 
in  aspect  angle  than  the  horizontal  case. 


When  the  angle  of  incidence  from  end-on  ia  30*,  Figures 
49-52,  the  vertically  polarized  cross  sections  are  in  close 
agreement  through  10*  of  bistatic  angle  as  before.  The 
lobes  of  the  horizontally  polarized  cross  sections  have 
about  a  1-2  dB  difference  to  about  10*  of  bistatic  angle, 
however,  the  nulls  are  significantly  different  at  5* 
bistatio  angle. 

Another  observation  can  be  made  regarding  Figures  47- 
65:  as  the  angle  of  incidence  approaches  broadside  the 

bistatic  angle  for  which  the  bistatic  and  equivalent 
monostatic  RCS  agree  becomes  larger.  This  concurs  with  the 
data  obtained  for  the  square  flat  plate.  For  both  vertical 
and  horizontal  polarizations  near  broadside,  the 
contribution  from  scattering  center  S3  is  negligible  while 
S|  and  S2  have  similar  diffraction  coefficients  and  phase. 
This  results  in  low  sensitivity  to  changes  in  aspect  angle 
and  good  agreement  between  bi static  and  equivalent 


monostatic  RCS. 
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Figure  47.  Cylinder  (•  =  15*;  Vertical  Polarization): 

Bi static  RCS  and  Equivalent  MonoBtatic  RCS 
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Figure  48.  Cylinder  (•  =  15*;  Vertioal  Polarization): 

Difference  in  RCS  (Biatatic  -  Monoetatic) 
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Figure  49.  Cylinder  (e=15#;  Horizontal  Polarization): 

Bi8tatic  RCS  and  Equivalent  Monoatatic  RCS 


Pigure  50.  Cylinder  (•=15*;  Horizontal  Polarization): 

Difference  in  RCS  (Bietatic  -  Monos tat io) 
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Figure  51.  Cylinder  (•  =  30*;  Vertical  Polarization): 

Biatatio  RCS  and  Equivalent  Monoatatic  RCS 


Figure  52.  Cylinder  (•  =  30*;  Vertical  Polarization): 

Difference  in  RCS  (Biatatio  -  Monoatatic) 
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Figure  53.  Cylinder  ( •=  30*;  Horizontal  Polarization): 

Bi static  RCS  and  Equivalent  Monostatio  RCS 


Difference  in  RC3  (Bistatio  -  Monoatatio) 
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Figure  55.  Cylinder  (•  =  45*;  Vertical  Polarisation): 

Biatatic  RCS  and  Equivalent  Monoatatic  RCS 


I 

Figure  56.  Cylinder  (•  =  45*;  Vertical  Polarisation): 
|  Difference  in  RC8  (Biatatio  -  Monos tat io) 
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Figure  57.  Cylinder  (e=45 * ;Horisontal  Polarisation): 

Biatatic  RCS  and  Equivalent  Nonostatio  RCS 


I 

I 

I 


Figure  58. 


Cylinder  (•=45*;  Horisontal  Polarisation) 
Difference  in  RCS  (Bistatio  -  Monostatic) 
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Fifure  59.  Cylinder  (•  s  60*;  Vertical  Polarisation): 

Bistatic  RCS  and  Equivalent  Monostatic  RCS 


Figure  61.  Cylinder  (e=60* ;Horisontal  Polarisation) : 

Bietatio  RCS  end  Equivelent  Monoatatio  RCS 
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Figure  63.  Cylinder  (•  =  7i5*;  Vertical  Polarisation): 

Biatatio  PCS  and  Equivalent  Monostatio  RCS 


Figure  64.  Cylinder  (•  s  75*;  Vertical  Polarisation) 
Difference  in  RCS  (Bietatic  -  Monostatio) 


Figure  65.  Cylinder  (•=76*;  Horisontal  Polarisation): 

Bi static  RCS  and  BQuivalent  Monos tat ic  RCS 
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Figure  66. 


Cylinder  (»iT5*;  Horisontal  Polarisation) 
Differenoe  in  RCS  (Bistatio  -  Monos tatio) 
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11  can  be  seen  in  Figures  67  and  68  that  the  sale 
problems  occurring  at.  end-on  also  occurred  at  broadside. 
The  only  difference  between  the  two  is  that  broadside 
biatatic  RCS  uses  eq.  (33)  while  end-on  uses  eq.  (38). 
There  is  little  difference,  though,  in  these  equations  and 
the  overall  effect  is  the  saae. 


I 
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Figure  67.  Cylinder  (•  *  90*;  Vertical  Polarisation) : 

Bistetio  BC8  end  Bquivalent  Monostotic  RCS 
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Figure  68. 


Cylinder  (6*90*;  Bo ri son tel  Polerisetion): 
Bistetio  BC8  end  Bquivelent  Monostetio  RCS 

TO 


IV. 


Conclus i on 


The  purpose  of  this  study  was  t.o  investigate  the 

applicability  of  Robert  Kell’s  method  under  varying 

■■  'll  ' 

conditions.  The  mathemathical  formulation  of  Kell’s  . 
relationship  was  first  examined.  This  was  done, by 
presenting  the  proof  for  the  monostatic-bistatic  equi Valence 
theorem  (MBFT),  the  foundation  of  Ke]  j*  s.  method",  and 
research  supporting  the  MBET.  It  was  brought  out  that,  the 
MBET  is  a  physical  optics  approximation  and,  therefore,  it. 

is  limited  by  physical  optics  requirements.  The 

■  '  ./ 

mathematical  formulation  of  Kell's  method  was  then 

presented.  With  this  formulation  it.  was  shown  Kell’s  method 
does  not  rely  on  approximations  of  the  scattered  fields  but 
upon  the  scattering  renters  of  the  target.  The  accuracy  of 
Kell’s  method  is  dependent,  upon  the  behavior  of  the 
individual  scattering  centers  relative  to  each  other  as  the 
hi  static  angle  changes. 

Kell’s  method  was  theoretically  applied  to  three  shapes 
of  varying  complexity  to  determine  some  of  its  limitations. 
The  sphere  was  used  to  show  how  well  Kell’s  method  predicted 
bi static  RCS  for  electrically  email  targets  as  compared  to 
large  targets.  Also,  the  effects  of  creeping  waves  were 
examined.  A  square  flat  plate  was  analyzed  to  investigate 
the  effects  angle  of  incidence  have  upon  the  accuracy  of 
Kell’s  method.  Finally,  a  right  circular  cylinder  was 
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studied  to  determine  what  success  Kell's  method  would  have 
predicting  the  bistatic  RCS  of  a  target  containing  surface 
discontinuities . 

Cone] us i ons 

This  investigation  was  conducted  to  determine  the 
applicability  of  Kell's  method  as  certain  parameters  were 
changed.  These  parameters  are  bistatic  angle,  electrical 
size  of  the  target,  angle  of  incidence,  and  surface 
continuity.  The  following  conclusions  were  reached. 

Ri  st.at  i  o  Angl  e .  Kell  claims  his  relationship  works 
best  when  the  scattering  centers  are  insensitive  to  changes 
in  bi static  angle.  The  scattering  centers  of  an 
electrically  large  sphere  are  insensitive  to  changes  in 
bi static  angle  relative  to  one  another  at  high  frequencies. 
As  a  result  Kell's  relationship  proved  to  be  quite  accurate 
for  spheres  at  high  frequencies  beyond  bistatic  angles  of 
80*.  The  scattering  center  of  the  flat  plate  and  circular 
cylinder  are  least  sensitive  to  changes  in  bistatic  angle  at 
high  frequencies  near  broadside  incidence.  Under  these 
condition,  Kell’s  relationship  was  accurate  (less  than  2-dB 
difference  between  bistatic  and  equivalent  monostatic  RCS) 
to  20*  of  bistatic  angle  for  the  square  flat  plate  and  30* 
for  the  cylinder. 

Electrical  size.  The  sphere  was  the  target  shape  used 
to  predict  the  accuracy  of  Kell’s  method  when  the  electrical 
size  of  the  target  is  varied.  It  can  be  seen  in  the  data 
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that  good  accuracy,  less  than  1 -dB  difference  for  both 
horizontal  and  vertical  polarization  did  not  occur  until  the 
sphere  became  large  (ka  =  10).  Much  of  the  inaccuracy  for 
smaller  spheres  can  be  attributed  to  the  inability  of  Kell’s 
method  to  account  for  creeping  waves  which  have  small  effect 
on  the  RCS  of  electrically  large  spheres. 

Angle  of  incidence.  It  can  be  seen  in  both  the  flat 
plate  and  cylinder  data  that  the  bistatic  angles  for  which 
the  bi static  and  equivalent  monostatic  RCS  agree  decreases 
as  the  angle  of  incidence  increases  from  broadside.  For  the 
flat  plate  some  of  the  degradation  in  agreement  between  the 
bi static  and  equivalent  monostatic  RCS  is  attributable  to 
the  inability  of  physical  optics  to  accurately  model  RCS  as 
the  angle  of  incidence  from  broadside  becomes  larger  than 
45*.  In  the  case  of  the  cylinder  using  horizontal 
polarization,  the  difference  between  bistatic  and  monostatic 
RCS  is  due  to  the  interaction  of  the  diffracted  fields  from 
the  scattering  centers.  As  the  aspect  angle  approaches  end- 
on  to  the  cylinder  the  diffracted  field  becomes  more 
sensitive  to  changes  in  aspect  angle  and  Kell’s  relationship 
loses  accuracy. 

Surface  continuity.  The  agreement  between  bistatic  and 
monostat.ic  RCS  becomes  polarization  dependent  when  surface 
discontinuities  in  the  form  of  90*  wedges  are  present.  The 
relative  phase  between  scattering  centers  was  sensitive  to 
changes  in  aspect  angle  when  horizontal  polarization  was 
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utilized.  This  resulted  in  poor  agreement,  between  the 
bi static  and  equivalent  monostatic  RCS .  Because  one 
scattering  center  dominated  the  RCS  when  vertical 
polarization  was  employed,  the  relative  phase  between 
scattering  centers  was  insensitive  to  changes  in  aspect 
angle  giving  better  agreement  between  the  bistatic  and 
equivalent  monostatic  cases.  Overall,  Kell’s  relationship 
produced  better  agreement  between  bi static  and  monostatic 
RCS  when  the  surface  was  electrically  large  and  free  of  any 
surface  discontinuities  such  as  a  sphere’s  than  when 
significant  discontinuities  such  as  90*  wedges  are  present. 
Recommendations  for  Further  Study 

Other  avenues  of  study  to  determine  the  limits  of 
Kell’s  method  include: 

1.  A  comparison  of  measured  bistatic  RCS  data  with 
measured  monostatic  RCS  data  at  appropriately  adjusted 
frequencies  is  a  required  step  in  the  process  of  verifying 
Kell’s  relationship.  Initially,  electrically  large,  simple, 
and  continuous  shapes  should  be  measured.  These  shapes 
reduce  higher  order  effects  and  their  RCS  is  not  affected  by 
diffraction  from  discontinuities.  The  complexity  of  the 
shape  can  then  be  increased  to  include  edges  such  as  on 
plates,  discontinuous  surfaces  such  as  the  mating  of  a 
spheroid  and  cylinder,  and  multiple  reflections  like  those 
found  in  a  hollow  cylinder  or  corner  reflector.  Also, 
various  materials  and  polarizations  could  be  measured  using 
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the  same  approach. 

2.  Using  the  approach  this  study  followed,  either  more 
rigorous  computer  modeling  of  the  simple  shapes  used  in  this 
investigation  or  computer  modeling  of  shapes  resembling 
aircraft  such  as  cylinders  with  hemispherical  or  spheroidal 
end-caps  could  be  examined.  The  RATSCAT  has  some  bistat.ic 
RTS  measurement  data  of  such  shapes.  Comparing  their  data 
with  data  derived  using  Kell's  relationship  using  computer 
modeling  would  give  a  good  indication  of  the  accuracy  of 
Kell's  relationship. 

3.  Study  the  corresponding  behaviors  of  the  monostatic 
and  bistat.ic  RCS  keeping  the  bisector  of  the  bi  static  angle 
oriented  at  a  constant  aspect  angle  to  the  target  as  the 

bi static  angle  increases.  This  approach  could  be  used  to 
study  the  accuracy  of  Kell’s  relationship  when  viewing  a 
particular  surface  discontinuity  or  to  present  the  data 
provided  in  a  study  such  as  this  one  in  a  different  light. 
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Appendix  A 


The  following  derivation  presents  the  mathematical 
formulation  Kell  employed  in  an  effort  to  relate  bistatio  RC8 
to  monos tat io  RCS . 

Radar  cross  section  can  be  expressed  in  terms  of  the 
Poynting  vector  (Kerr,  1951:33): 


2|s's 

o  =  4kRz  — 

Is1 

o  =  radar  cross  section 

Ss  =  back  scattered  Poynting  vector 

S1  =  incident  Poynting  vector 


(34) 


For  I  E  |  =  |  H  |  ( u0  /€0)^1/Z)  and  taking  the  time  average 
value  of  S,  eq.  (34)  becomes 

S  =  ( 1/2 ) |E | |H |  (35) 

av 


and,  therefore, 


H8  2 ] / 2 


o  =  4*R 


2  -  *  ^o'  o 


[(Mo/Co)^/2)  Ih1  1 2 3 / 2 


(36) 


Reducing  and  assuming  far  field  gives  the  starting 
point  for  Kell’s  formulation: 


o 


4n  lim  RQ 
R0-*« 


(37) 
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where  Hp  =  H8  =  reradiated  magnetic  field 
R0  =  distance  to  observer 

Solve  for  Hp  using  the  Stratton-Chu  equation  (Stratton, 
1941:466): 


Ho  = 


1 

Ttx 


J  x  VV  +  “^mVvldV 


I  ( i w€M  + 

J\f 

-  4ii  f  t(fi  x  "8>  x  V*>  +  (fi  X  Hs)Vv 


-iu>€(fi  x  Er)v»]  da 


(38) 


M  =  magnetic  current  density 
J  =  electric  current  density 
m  =  magnetic  charge 
Hs  =  magnetic  vector  on  target  surface 
Es  =  electric  vector  on  target  surface 

J 

i>  -  free  space  Green's  function  (e  °/r0; 

rQ  =  distance  from  elemental  area  to  observer) 

If  we  assume  there  are  no  sources  within  the  target  volume 

then  Stratton’s  equation  can  be  reduced  to  give  the  right 

hand  side  of  Kell’s  equation  (2)  given  below  (Kell, 

1965:983) : 

N"Hr  =  -  4“  J  t(rt  x  Ha)  x  7eikro  +  (fi  x  )  V eikro 
s  ro  ro 

-iw€ ( A  x  Es)(ekro)/r0]  da  (39) 


Hp  =  magnetic  vector  amplitude  of  reradiated  field 

C  =  phase  of  reradiated  field  relative  to  chosen 
reference 
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da  s  element  of  area  on  surface  of  target 


ft  =  unit  vector  normal  to  target  surface  and 
pointing  inward 

Figure  69  relates  the  target  to  the  incident  and  scattered 
fields . 


Figure  69.  Biatatic  Coordinate  System  with  Antennas 
in  the  x-z  Plane  (Kell,  1965:984) 

Kell  gives  the  relation  between  rj,  rQ,  Rj,  RQ,  s  and 

B  as 

(r,+r0)  =  2zCos ( i/2 )  4  2(R|4R0)  (40) 

This  is  incorrect.  The  correct  relationship  is 

(rj+r0)  s  2sCos  ( i/2 )  4  (Rj4R0)  as  Rj  and  R0  (41) 
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Kell  normalizes  the  surface  fields  by  dividing  by  the 
absolute  magnetude  of  the  incident  field.  He  then  accounts 


for  illumination  phase  delay  at  each  elemental  area  by 
dividing  by  ei(kri+*)  ,  This  is  summarized  in  eqs  (42)  and 


(43)  . 

For 

l*sl 

-i  (krj+.l  ,  -  and  l^-Ukr,**) 

IejI 

8  Iff,  1 

|ES| 

a  es  |E|  |  exp(  i  (k^t*)  ) 

IhJ  =  hs  |Hj  I  exp{  i  (krj+*) ) 

where 


(42) 

(43) 


eg  =  normalized  electric  vector 

hs  =  normalized  magnetic  vector 

*  =  phase  difference  between  the  incident  field 
(either  electric  or  magnetic)  and  the  local 
field  at  da 


Using  far  field  approximations 

V«kr°  -  Ul  („ikro)  p 
(44)  "  '*  ' 

ro  ro 

And  recognizing 

I  I’ll  =  I  Hi  |  ( P0/  €c)  (4/2)  (45) 

Equations  (42)-(4S)  can  be  substituted  into  (39)  to 
give 

V«R  =  -4“  X  [  < ft  X  irs|Hi|ei(kri+#)  )  x  fQ(  ikeikro  ) /Tq 
+  <fi.h8|Hj|e  i(kri+*))f0(ikeikro)/r0 

-  i w(m0/€0)  1/2(A  x  e8|Hi|e  1  (kri+*>)  (e  ikr0)/ro]da  (46) 
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Rearranging  the  above  equation  gives 

fo"»R  s  x  h8ei(kri4#)>  x  fQ(  ikeikro  > /rQ 

+  (ft.h8  «i(kri+#)  )f0(ikeikro)/r0 

-  i«(M0/€0)  1/2(ft  x  eae  i(kri+4)  )(eikro)/r0]da  (47) 

Knowing  that  w(m0/«0)*^  *  k  and  making  the  far  field 
approximation  rQ  ~  RQ  ,  eq.  (47)  can  be  written 

fo'HR  =  *  h.«lk,ri+ro)  >  *  *<,«**  > 

*  i(n.h8eik(ri+ro>  )(e*)tQ 

-  i  (ft  x  e8eik(ri+ro))ei#  ]  da  (48) 

Using  the  relationship  given  in  eq.  (41),  the  phase  term 
in  eq.  (48)  can  be  rewritten  as 

eik(rj+r0)  _  ei2kzCos(8/2)  ^MR^Rq)  <49) 


x  k/  2 

By  writing  the  imaginary  expression  i  as  e  and  using  eq. 

(49),  eq.  (48)  becomes 


Hj 


*oHR  = 


2R0x 


ik(Ri+R0)  ix/2 


(ft  x  h8)  x  fr( 


A** 


+  (ft.h8)(ei#)f0  -  (ft  x  e8)ei#) 
x  eik2sCo.(l/2)  da 


(50) 


Any  specific  geometry  will  be  dependent  on  •  and  z, 
therefore,  eq.  (50)  can  be  rewritten  using  p(«,s),  a 


80 


conversion  factor  relating  the  surface  element,  da,  to  the 
differential  coordinate  product  dodz  to  yield  Kell’s  eq.  (6) 
with  one  difference,  the  ej  term: 


*o»R  = 


ik(Ri+R0)  i*/2 


2RnX 


J J  [  (  ft  x  h  8 


)  X  fn(J*  ) 


+  (fi.hs)(eU)fG  -  (fix  es)e  **) 

.  ik2zCos(0/2) 
x  pi®,*)  e  '  dedz 


(51  ) 


Using  eq.  (51)  in  (37)  yields 


_*  pi^(VRo>  e 
x2 


ia/2  jj 


t(fi  x  hB)  x  fiD(e  M 


+  ( n  *hs)  ( e  1  *)  fQ  -  (fix  e8)e  l' 


p(e,z)  eiR2^os(g/2)  dedz 


(52) 


By  allowing  the  magnitude  squared  of  e*k^Ri+Ro)  e*n^2 
equal  l  and  integrating  eq.  (52)  with  respect  to  6,  the  RCS 


is  given  by 


■f  |JlU)  .l*k»Co.(»/*)dE|  2 


(53) 


where 


;)  =  J[(fi  x  h8)  x  fG(ei#)  +  (fi.h8)(ei#|f0 

-  (ft  x  e8)ei#]  p(0,z>  eik2zCo8(fl/2)  do 


(54) 
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The  following  comments  can  be  made  regarding  eq.  (53): 

a)  It  has  the  form  of  physical  optics. 

b)  The  difference  between  eq.  (53)  and  physical 
optics  lies  in  the  fact  that  eq.  (53)  is  exact.  However, 
the  T(z)  term  in  eq.  (53)  is  not  precisely  known. 

c)  T ( s )  is  a  composite  of  surface  geometry, 
illuminating  and  observing  ray  geometry,  and  surface  wave 
propagation  effects. 

d)  The  analytic  continuity  of  T(s)  significantly 
effects  o  since  eq,  (53)  may  be  subdivided  into  a  sum  of 
integrals,  each  subintegral  taken  over  a  range  of  *  wi thing 
which  its  integrand  is  continuous. 

e)  Approximations  of  I(z)  yield  contributions 
dependent  only  on  the  integral’s  endpoints. 

a 

f)  The  endpoints  are  the  scattering  centers. 

From  this  logic,  Kell  proposes  that  RCS  can  be 

computed  from  discrete  scattering  centers  as  shown  below: 

o  -  |E  (om)l/2eipm  |  2  (55) 

m=  1 

M  =  number  of  scattering  centers 

om  =  RCS  of  m**1  scattering  center 

=  phase  of  field  scattered  by  m**1  center 
relative  to  the  phase  of  the  first 
scattering  center 

According  to  Kell,  eq.  (55)  is  used  to  determine  the 
relation  between  monostatic  and  bistatic  RCS  for  a  given 
aspect  angle  a.  Expanding  the  phase  term  in  eq.  (55) 
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<  56 ) 


*m  =  2k7.mCo«<*/2)  ♦ 


where  em  = 
*m  s 


f  ubst  i  i  ng 


d i stance 

residual 
etc,  ) 

eq .  ( 56  ' 


between  l8*1  and  m1^  scattering  center 
phase  contributions  (creeping  waves, 

into  ( 55 ) : 


£  (om)1/2  ei2k2mCos<P/2> 


(57) 


Kell  uses  eq.  (57)  to  make  a  statement  for  monostat.io- 
bi static  equivalence  and  the  conditions  required  for  its 
occurrence : 

Tf  for  a  chosen  aspect  angle,  the  following 
conditions  hold: 

1)  The  RCS  may  be  written  as  a  squared  sum  of  fields 
from  discrete  scattering  centers,  and 

2)  The  amplitude  •  position  zm,  and  residual 

phase  are  insensitive  to  the  bistatic  angle  B  over 
the  range  of  B  considered,  for  those  centers  which  are 
significant  members  in  this  sum; 

it  then  follows  that  the  bistatic  cross  section  of 
aspect  angle  ®  and  bistatic  angle  B  is  equal  to  the 
monostatic  cross  section  measured  on  the  bisector  at  a 
frequency  lower  by  the  factor  Cos (0/2)  [Kel 1 , 1 965 : 987 ] . 
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Annend 1  x  B 


This  nppendix  present *  the  computer  programs  used  to 
calculate  the  bi static  and  equivalent  monostat.ic  RCS  of  a 
sphere.  These  programs  were  developed  by  Ohio  State 
University  in  197ft.  They  were  then  modified  by  AFWAL  for 
use  on  the  indoor  radar  range  at  Wright-Patterson  AFB.  The 
first  programi  BTSPH.FOR,  computes  the  bistatic  RCS. 
MONOSPH.FOR  is  a  modification  to  BTSPH.FOR.  Tt.  calculates 
monostatic  RCS  at  a  frequency  adjusted  by  a  factor  of 


Sec(0/2)  (0  =  bistatic  angle)  to  give  a  one-to-one 

correspondence  between  bi static  RCS  and  monostatic  RCS. 

C  ================<<<  BTSPH.FOR  >>>==================== 

C 

C  THIS  PROGRAM  GENERATES  SPHERE  CALIBRATION  DATA  FOR 
C  THE  SWEPT  FREQUENCY  SYSTEM 
C 

C  OUTPUT— MAGNITUDE:  10LOG(RCS  )-40 ;  RCS  IN  SQUARE  CM 

C  PHASE:  RADIANS;  FROM  SPHERE  CENTER 

C 

C  MODIFICATIONS:  THIS  PROGRAM  WAS  MODIFIED  FOR  USE  ON  A 
C  PDP- 11/23.  OUTPUT  IN  NOW  STORED  IN  ONE 

C  2-DIMENSIONAL  ARRAY  YM.  THE  PHASE  IS 

C  GIVEN  IN  RADIANS.  OUTPUT  GOES  TO  AN 

C  UNFORMATTED  FILE. 


C 

COMMON  BUFF,NDIM, ANST , AINC 
COMPLEX  ETH.EPH 
character*20output  file 
REAL  KA 

DIMENSION  AM ( 2000) , PH ( 2000) 
INTEGER* 2  INFILE<15) 

BYTE  BUFF (35000) 

LOGICAL  VHP.LBK 
'DATA  PI/3.141593/ 

DATA  LT,LF/-1 ,0/ 

CONST: 2 . *PI *  1 . E9/300 . E6 
RTD= 180 . /PI 
C 
C 

C  ====DKFINE  RANGE  GEOMETRY 
C  SPHERICAL  SCATTERING 
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WRITE ( S , * )  ’INPUT  T  FOR  VER  POL  AND  F  FOR  HOR  POL’ 
ACCEPT  *,  VHP 

WRITE (6,*)  ’INPUT  ANTENNA  START  SEPARATION  IN  DEGREES’ 

ACCEPT  *,  ASD 

asdpsasd 

write (6,*)  ’ TNPUT  ANTENNA  STOP  SEPARATION  IN  DEGREES’ 

accept,  *,  asdl 

THRzlfiO .-ASD 

THBp=thb 

PHR=-90 . 

IF (VHP)  PHR= 1 80 . 

write (6,*)  'ANTENNA  SEPARATION  INCREMENT  >=  1  DEG’ 
accept.  *  ,dela 
nba= ( asd 1 -asd ) /del at  1 .01 
C 

WRTTE ( 6 , * )  ’INPUT  SPHERE  DIAMETER  IN  INCHES’ 

ACCEPT  *,SDT 

WRTTE (6,*)  ’INPUT  MIN  FREQ.  (GHZ):  ’ 

ACCEPT  *,  FMIN 
ANSTsFMIN* 1 000  . 

WRITE ( 6 , * )  ’INPUT  MAX  FREQ.  (GHZ):  ’ 

ACCEPT  *,  FMAX 

WRITE ( 6 , * )  ’FREQUENCY  INCREMENT  (GHZ):  ’ 

ACCEPT  *,  DELF 

write(6,t)  'Type  in  output  filename’ 
read  (5,  15)  out.put_file 
15  format (a20) 

NF= ( FMAX-FMTN ) /DELF+ 1 . 1 
NDTM=NF 
C 

SRCM=SDI*2 . 54/2. 

SRM=SRCM/100. 

AREACM=4 . *PI*SRCM**2 
AINC=DELF* 1000. 

CKA= CONST *SRM 
FREQ=FMIN-DELF 

=r==START  OF  LOOP==== 


FUNCTION  CATAN2 ( Z ) 
COMPLEX  Z 
RZsREAL ( Z ) 
FIZ=AIMAG(Z) 

CATAN2= ATAN2 ( FIZ , RZ ) 

RETURN 

END 
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SUBROUTINE  OPS ( VHP , KA , EV , ETH , EPH ) 

COMPLEX  ETH , EPH , ETHT , EPHT , ETHI , EPH T 

REAL  KA 

LOGICAL  VHP 

IF( VHP)  THEN 

THT=180. 

PHT=180. 

CALL  FIELD ( 1 . , KA , THT , PHT , ETHT , EPHT , IER ) 

ETH=ETHT 

THT=180.-2*EV 

PHT  =  1 80 . 

CALL  FT ELD ( 1 . , KA , TH I , PHI , ETH I , EPHI , TER) 

ETH=ETH+ETHI 

ELSE 

THT= 180. 

PHT  =  -90 . 

CALL  FIELD  (  T  .  ,  KA  ,  THT  ,  PHT  ,  ETHT  ,  EPHT  ,  IER  ) 

EPH=EPHT 

THT  •=180.-2.  *EV 

PHT  =90 . 

CALL  FIELD ( 1 . , KA , THT , PHI , ETHT , EPHT , TFR ) 

EPH = EPH + EPH I 

ENDTF 

RETURN 

END 


B I STATIC -BACKS CATTERED  FIELD  OF  A  SPHERE 

SUBROUTINE  FIELD ( TIMCON , KA , THE , PHI , GTHE , GPHI , IER ) 

TIMCONsTIME  CONVENTION=+ 1 . 0  FOR  HARRINGTON  ( - JWT ) 

-1.0  FOR  STRATTON  ( + JWT ) 

COMPLEX  J, GTHE, GPHI, SGTHE,SGPHI ,BN,CN 
DIMENSION  SJ( 1 50) ,SY( 150) ,DSJ( 150) ,DSY( 150) ,DP( 150) 
DOUBLE  PRECISION  TPCERR , ANGDIF , P ( 151 ) , DEF 
DOUBLE  PRECISION  U , DCOS , DSTN , V 
REAL  KA 
C 


C=====DEFTNE  CONSTANTS 
C  J =SQRT ( - 1 ) 

C  TPCERR = TOTAL  X  CHANGE  ERROR  ALLOWABLE  BETWEEN 

C  SUCCESSIVE  ITERATIONS  TO  DEVINE  CONVERGENCE 

C  TPCMAGsTOTAL  X  CHANGE  IN  MAGNITUDE  ACTUALLY 

C  OCCURRING  BETWEEN  SUCCESSIVE  LOOPS. 

C 

DATA  PI, J/3. 141593, (0. , 1 . )/ 

DATA  TPCERR, ANGDIF/1 .OD-20, 1 .OD-2/ 

C 

C===== INITIALIZE  VARIABLES 
C 
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GTHF,  =  CMPLX  (  0  .  ,0.  ) 

GPHT  =CMPI.X  (  0  .  ,0.  ) 

SGTHEsCMPLX ( 0 . ,0. ) 

SGPHT  sCMPLX ( 0 . ,0. ) 

C 

MO= 150 

CALL  SPHEBES(KA,SJ,SY,DSJ ,DSY,MO,MAX) 

C 

Cs=s=s ELIMINATE  ZERO  ORDER  TERMS  FROM  THE  ARRAYS 

C 

MAX=MAX- 1 
DO  5  Ls 1, MAX 
LL-I.+  1 

S  J  ( I. )  =S  J  ( I.L  ) 

SY(L)=SY(LL) 

DS  J  ( I, )  =DSJ  {  LL ) 

DSY ( L ) sDSY ( LL ) 
fi  CONTINUE 

C 

CHKTHE=ARS (180. -THE ) 

IF  ( CHKTHE  .LT.  ANGDIF )  GOTO  20 
IF  (ABS(THE)  .I.T.  ANGDIF)  GOTO  15 
U=DCOS ( DBLE ( THE#PI /I 80 . ) ) 

Vs DSTN ( DBLE ( THE* PI / 1 80 . ) ) 

MAXNsMAX+1 
M=  1 

DEF= - 1 .DO 

CALL  POLY2 { DEF ,  U ,  M , MAXN , P ) 

DO  10  Nsl , MAX 

DP ( N  )  s  ( FLOAT (N+1)*U*P(N) -FLOAT ( N-M+ 1 ) *P(N+1 ) 
&*V/1 . -U*  *  2 ) 

P(N)=P(N)/V 
10  CONTINUE 

GOTO40 

15  DO  16  Nsl, MAX 
FN=FLOAT(N) 

DP(N)  =  .5*FN*(FN+1  . ) 

P(N)=(— 1 . )*DP(N) 

16  CONTINUE 
GOTO  40 

20  DO  30  Nsl , MAX 

FNsFLOAT(N) 

P(N)s0.5*FN*(FN+l  )*(-l  .  )MN 
DP(N)=P(N) 

30  CONTINUE 

C 

40  DO  50  Nsl , MAX 

AN=  (  -1  .  )  *  FLOAT  (  2  *  N+ 1  )  /  FLOAT  ( NM2  +  N ) 

BNs AN* ( S J ( N ) +KA*SDJ ( N ) / ( S J ( N ) +KA*DS J ( N ) - J* 

& ( SY ( N ) +KA*DSY ( N ) ) ) 

CNsAN*SJ(N)/(SJ(N)-J*SY(N) ) 

GTHEsGTHE+BN*DP(N)-CN*P(N) 

GPHIsGPHI+BN*P(N)-CN*DP(N) 
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nan 


AA=CARS ( GTHE-SGTHE ) 

BBsCARS ( GTHE ) 

CC=CABS ( GPHI-SGPHT ) 

DDsuABS  (GPHI  ) 

IF  (BB  .NE.  0. )  GOTO  44 
AA  =  0  . 

RR=]  . 

44  IF  (DD  .NE.  0. )  GOTO  45 

CC  =  0.  -v 

DD=  1  . 

-15  TPCMAG=AA/BR+CC/DD 

TER  =  N 

TF  (TPCMAG  .LF.  TPCERR )  GOTO  60 
SGTHE=GTHF 
SGPHT  =GPHT 
50  CONTINUE 

TER  =  0 

r 

60  GTHE=TI MCON*GTHE* J*COS { PHI #PI / 1 80 . ) /KA 

GPHT  =TIMCON*BPHT * J*SIN ( PHI  API/ 1 80 . ) /KA 
RETURN 
END 


SUBROUTINE  SPHEBES ( X , B J . BY , BP , YP , I DM , MAX ) 

C 

C  X= ARGUMENT 

C  B»T  =  SPHERICAL  BESSEL  FUNCTION  ARRAY 

C  RY=SPHERICAL  NEUMAN  FUNCTION  ARRAY 

C  BP=PRTMED  BESSEL  FUNCTION  ARRAY 

C  YP=PRIMED  NEUMAN  FUNCTION  ARRAY 

C  IDM=MAX  NUMBER  OF  ORDERS  TO  BE  COMPUTED 

C  MAXsMAX  NUMBER  OF  ORDERS  ACTUALLY  COMPUTED 

C 

C  COMPUTATION  IS  STOPPED  WHEN  THE  VALUE  OF  A  NEUMAN 

C  FUNCTION  MAGNITUDE  IS  GREATER  THAN  FMAX  WHICH  IS 

C  DEFINED  INTERNALLY  TO  THE  PROGRAM. 

C 

DIMENSION  B J ( 1 ) , BY ( 1 ) , BP ( 1 ) » YP ( 1 ) 

FMAXs 1 . E35 
FMIN=1 .E-38 
SX=SIN ( X ) /X 
CX=COS(X)/X 
FF1=SX 
FF2=SX/X-CX 
VY1=-CX 
Yf2=- ( CX/X+SX ) 

BY( 1 )=YY1 
BY ( 2 ) =YY2 
BYB=YY1 
BYC=YY2 
1  =  2 
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20  BYA=BYB 
BYB=BYC 

BYCs ( 2 . *  I- 1 . )*BYB/X-BYA 
L=  1  + 1 

IF  (L  .IE.  IDM)  BY ( L ) =BYC 
AY1=ABS(BYC) 

1  =  1  +  1 

IF  ( AY1  .LT.  FMAX )  GOTO  20 

MAX=I -  I 

BJBsO. 

BJAsFMIN 
I  =MAX 
50  1=1-1 

BJC=BJB 

BJB=BJA 

BJA= ( 2 . *  1  +  3 ) *RJB/X-BJC 
L=T  +  1 

IF  (L  .LE.  IDM)  BJ(L)=BJA 
IF  (I  .GT.  0)  GOTO  50 
ALF=BJ(2)/FF2 

IF  (ABS(FFl)  .GT.  ABS(  FF2  )  )  ALF=BJ ( 1 ) /FF1 

ALF=1./ALF 

K=MAX 

IF  (K  .GT.  IDM)  K=IDM 
DO  60  1  =  1, K 
60  BJ(I)=BJ(I)*ALF 

BP( 1 ) =-BJ ( 2 ) 

YP( 1 )=-BY(2) 

DO  80  1=2, K 
I  ri  =  I  —  l 
FAC=I/X 

BP( I ) =BJ ( IM ) -FAC*BJ ( I ) 

80  YP( I ) =BY ( IM ) -FAC3BV ( I ) 

RETURN 

END 

C 

C============================== ==--================ 

c 

C  ASSOCIATED  LEGENDRE  POLYNOMIAL  SUBROUTINE 

C 

C  FOR  HARRINGTON'S  DEFINITION  OF  THE  ASSOCIATED 
C  LEGENDRE  PLOYNOMIAL  DEF=-1. 

C 

C  FOR  THE  OTHER  DEFINITION  OF  THE  ASSOCIATED 

C  LEGENDRE  POLYNOMIAL  DEF=1. 

C 

SUBROUTINE  POLY 2  ( DEI  , ?. ,  M ,  MAXN ,  P ) 

DOUBLE  PRECISION  DEF , X , P ( 1 ) , SQ , DSQRT , DBLE 
DOUBLE  PRECISION  FL1 , FL2 , FL3 
SQ=DSQRT( 1 . 0D0-X**2 ) 

P( 1 )=1 .DO 

IF  (M  .EQ.  0)  GOTO  1 
DO  2  L=1,M 
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P(  1  ) =DEF*DBLE( FLOAT ( 2*L- i ) ) *SQ*P( 1 ) 

2  CONTINUE 

1  P( 2 ) sDBLE ( FLOAT ( 2*M+ 1 ) )*X*P(1) 

DO  3  K=3,MAXN 
I  =K- 1 
J = K  -  2 
N=M+K- 1 

FL1 sDBLE ( FLOAT ( N  +  N-l )  ) 
FL2=DLBE(FLOAT(M+N-l ) ) 

FL3 sDBLE ( FLOAT { N-M ) ) 

P(K)s(FLl*X*P(l )-FL2*P(J) )/FL3 

3  CONTINUE 
RETURN 
END 


FUNCTION  ACOS(X) 

DOUBLE  PRECISION  DSQRT 
DATA  PI/3.141593/ 

IF  (X  . GE .  1.)  GOTO  1 

IF  (X  .LE.  -1.)  GOTO  2 

ACOSs ( PI/2 . )-ATAN(X/DSQRT( 1 . D0-X4X) ) 

RETURN 

ACO^O. 

RETURN 

ACOSrPI 

RETURN 

END 


The  following  routine  is  the  difference  between 
.FOR  and  MONOSPH.FOR. 


===<<<  MONOS PH. FOR> >>==== 


BISTATIC  SCATTERING  FROM  A  SPHERE 

open ( uni t=10,'.<tatU8=  'NEW*  ,  filesoutput_file, errs  100) 
writeUO,*)  *sph.  Sep.  Angle  freq  Mono.' 

write (10,*)  ’size  (deg)  (GHz)  RCS * 

write(  10 ,  *  )  Min)1 
DO  45  jjsl.nba 

aad=asdp+ j j-1 
thb=180 . 

fref;^fain»l/(co8(  .  5*asd/rtd)  ) 

KAsCKAt FREQ 

CALL  FIELD ( 1 . , KA , THB , PHB , ETH , EPH , TER ) 
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RCS  FOR  HORIZONTAL  POLARIZATION 

IF(VHP)  THEN 

RCS = ARE ACM* CABS ( ETH ) *  *  2 
PH ( I ) =CATAN2 ( ETH ) *RTD 

RCS  FOR  VERTICAL  POLARIZATION 

ELSE 

RCS = ARE ACM* CABS ( EPH )  *  *2 
PH ( I ) =CATAN2 ( EPH ) *RTD 

ENDIF 

AM( I ) : 10 . * ALOG1 0 ( RCS ) -40 

wri 16(10,50)  sdi , asd , f req, AM( I ) 

5  continue 

cloae(uni 1=10) 
goto  101 

00  write(6,*)  ’Filename  already  exists!!!* 

0  FORMAT  ( f 5 . 2 , 2X , F10 . 4 , 2X,F1 0 *  4 , 2X ,F1 0 . 4 ) 
01  continue 
STOP 
END 
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APfltnflix.  C 


This  appendix  presents  the  derivation  of  the 

i 

mathematical  model  for  the  RCS  of  a  square  flat  plate  using 
physical  optica  approximations.  Following  the  derivation  is 
the  computer  programs  used  to  calculate  the  bistatic  and 
equivalent  monostatic  cross  sections  for  a  flat  plate. 

The  geometry  utilised  in  this  derivation  ia  given  in 
Figure  70. 


Figure  70.  Scattering  Geometry  for  a  Square  Flat  Plate 

The  incident  and  scattered  propagation  vectors  are 
constrained  to  lie  in  the  x-y  plane. 

The  incident  electric  field  can  be  expressed  as 
(Harrington,  1961:140) 

E1  s  E  e«J*{xCos#  4  ySina)  j58 


92 


r 


The  scattered  electric  field  can  be  expressed  as 
(Harrington,  1961:140) 

,  kK0a2e“Jkr  Sin[k(a/2 )  (Sin#  ♦  Sins)] 

Ej  m  Cos#  — -  — ■  ■' 

J2«r  k (a/2 ) (Sine*  Sine) 


From  eqs .  (35)  and  (36),  RCS  is  given  by 


=  4k  lim 

r—e*» 


Substituting  eqs.  (68)  and  (69)  into  eq.  (60)  yields 


(59) 


(60) 


O  S  4  K 


Cose 


Sin[k(a/2)(Sin»  »  Sinelll2 


: ( a/2 ) (S<n*  ♦  Sine) 


(61) 


From  Figure  (70),  it  can  be  seen  e  s  •♦8.  Substituting 
this  relationship  into  eq.  (61)  gives  the  expression  for 
bi static  RCS  in  terms  of  angle  of  incidence  and  bistatio 


ang] e : 


O  -  4k  I"—  Cos  ( e+ A )  Sin(k  ( a/2 )  ( Sine  ♦  Sin(»+g))n  2 
°  4*  L  x  Coa(#+fl>  k(a/2)(Sin«  ♦  Sin(e+8))  J  <62) 

In  accordance  with  Kell's  hypothesis,  the  aspect  angle  for 
the  equivalent  monostatic  case  is  equal  to  the  angle  of 
incidence  plus  half  the  bistatic  angle.  Applying  this 
ooncept  to  eq.  (61)  gives  -he  equivalent  aonostatio  RCS  as 


a  s  4k 


[jdCo.,..,/„s^§s| 


(63) 
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The  computer  programs  calculating  the  bistatic  and 
equivalent  monostatic  RCS  follow.  These  programs  are 
written  in  BASIC  and  run  on  a  Zenith  Z-100  computer  using  an 
8-bit.  CP/N  operating  system. 


1  >»»t*tMttt**tt*t**tBTSTATIC  RCStMtMtttttMMttMtM 

2  ’THIS  PROGRAM  CALCULATES  THE  BTSTATIC  RCS  OF  A  SQUARE 

3  *  FLAT  PLATE 

4  * 

8  ‘OUTPUT:  10*LOG(RCS);  RCS  TN  SQUARE  METERS 

6  ’ 

10  OPEN  "O" | 1 i "B: FPRCSB.DAT" 

1 1  * 

12  ‘DEFINE  VARIABLES 

13  ‘ L  =  LENGTH  OF  ONE  RIDE  IN  METERS 

14  *WL  =  WAVELENGTH  IN  METERS 

15  ‘THETA  =  ANGLE  OF  INCIDENCE 

16  ‘BETA  =  BISTATIC  ANGLE 

17  ’ 

20  L=0 . 3048 

30  PT=3. 141593 

40  DIM  RCSU81) 

44 

45  *T  =  ANGLES  OF  INCIDENCE 

46  ' 

50  FOR  1=0  TO  40  STEP  8 

60  A$  =  "ANGLE  OF  INCIDENCES " 

70  PR1NT#1,A$;I 

74  ‘ 

75  ‘N  =  BISTATIC  ANGLE 

76  ‘ 

80  FOR  N  =  0  TO  60 

90  THETA  =  PI *1/180 

100  BETA  =  NtPI/180 

110  WL=0.03 

140  J  =  L*PI *( SIN ( THETA ) 4SIN ( THETA+BETA ) )/WL 

150  IF  J  =  0  THEN  16*  ELSE  160 

160  RCS(N)=4*PI#(L*2*COS(THETA+B3TA)*SIN(J)/( J»WL) )*2 

161  GO  TO  170 

162  RCS(N)  =  4*PI* ( L* 2* COS (THETA+BETA) /WL ) *2 

170  NEXT  N 

190  FOR  K=0  TO  60 

210  RCS (K ) = 1 0*LOG (RCS (K ) )/2. 3025851 

220  PRINT*1 .USING"####.###  ";K,RCS(K) 

230  NEXT  K 

235  NEXT  I 

240  CLOSE# 1 

250  STOP 

260  *«*»****M#####*****t ********************************* 
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1  >tM**ttt*»t*t**EQUlVALBNT  NONOSTATIC  RCS*M********** 

2 

3  'THIS  PROGRAM  CALCULATES  THE  MONOSTATIC  RCS  OP  A  SQUARE 

4  'FLAT  PLATE  AS  A  FUNCTION  OF  FREQUENCY  IN  ACCORDANCE 

5  'WITH  KELL'S  HYPOTE8IS 

6  ’ 

7  'OUTPUT:  10*LOQ(RCS);  RCS  IN  SQUARE  METERS 

ft  » 

10  OPEN  "0" , 1 , "B: FPRCSF.DAT" 

11  ' 

12  'DEFINE  VARIABLES 

13  *L  a  LENGTH  OF  SIDE  IN  METERS 

14  *WL  =  WAVELENGTH  IN  METERS 

15  'THETA  =  MONOSTATIC  ASPECT  ANGLE 

16  'BETA  a  BISTATIC  ANGLE 

1?  » 

20  Ls0.3048 

30  PI=3. 141593 

40  DIM  RCS ( 181 ) 

45  DIM  FREQ(IOO) 

4?  ' 

48  'I  a  ANGLE  OF  INCIDENCE  FOR  BISTATIC  CASE 

49  ' 

50  FOR  1=0  TO  40  STEP  8 

60  AS  =  "ANGLE  OF  INCIDENCE:" 

70  PRINT#1 , AS ; I 

74  * 

75  *N  =  BISTATIC  ANGLE  IN  DEGREES 

76 

80  FOR  N  =  0  TO  60 

90  THETA  *  PI* ( I+N/2 )/180 

100  BETA  =  N*PI/180 

110  WL=0.03*COS( BETA/2) 

130  FREQ(N)  =  3B8/(WL*1E9 ) 

140  J  =  L*PI*2*SIN( THETA )/WL 

150  IF  J  =  0  THEN  162  ELSE  160 

160  RCS (N ) =4*PI f (L"2*COS ( THETA) *SIN ( J )/( J*WL ) ) *2 

161  GO  TO  170 

162  RCS(N)  =  4*PI*(L"2*COS(THETA)/WL)"2 

170  NEXT  N 

190  FOR  K=0  TO  60 

210  RCS(K)slO*LOG(RCS(K))/2. 3025851 

220  PRINT# 1 (USING"#### « ###  " ;FREQ(K) , RCS(K) 

230  NEXT  K 

235  NEXT  I 

240  CLOSE# 1 

250  STOP 

260  ****************************************************** 
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Appendix  J) 


This  appendix  presents  the  derivation  of  the  RCS  model 
for  a  right,  circular  cylinder  and  the  computer  programs  used 
to  calculate  its  histatic  and  equivalent  monostatic  RCS. 

The  derivation  of  the  RCS  is  taken  from  "Project  DISTRACT— 
Distributed  Radar  Concepts  and  Techniques, "  Appendix  A, 
Cornell  Aeronautical  Laboratory,  Inc. 

The  RCS  of  a  finite  cylinder  ia  dominated  by  diffracted 
returns.  For  this  reason,  geometrical  theory  of  diffraction 
(GTD)  is  used  to  compute  the  RCS,  A  general  diffraction 
coefficient,  is  written  ns 

|  TT  /  4 

D  =  2 - SiMtt/n)  [(Co8(i»/n)-Co*|«'-«*))‘1 

nsine'(2*k) 

?  (Cos( n/n)-Cos( ( •'♦o'+n)/n) )-1 ]  (64) 


where  o*  =  angle  of  incidence 

•*  =  angle  of  diffraction 

n  a  T/s  (T  s  exterior  wedge  angle) 

=  angle  between  the  incident  ray  and  the  positive 
tangent  to  the  edge  (*'  =  k/2  in  this  appendix), 

k  =  2*/X  (wave  number) 

Polarization  determines  the  choice  of  signs  used  in  eq.(64). 
Vertical  polarization  (E  parallel  to  the  edge  of  the  wedge) 
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uses  the  upper  sign.  Horizontal  polarization  uses  the 
lower.  Angles  ®'  and  a'  are  defined  graphically  in  Figure 
71.  Because  a  right  circular  cylinder  has  T  =  270*, 


n  -  3/2.  With  e'  =  */2,  eq.(64)  becomes 


D  =  |  el--4Sin  (2n/3)  t  {Cos  (  2n/3  ) -Cos (  e'-a ' )  ) 
3  (2nk)1/Z 


*  (Cos  (  2n/3  )  -Cos  (  2  ( ®  '  +  a'  +  n)  /?. )  )  ~ 1  ]  (65) 


INCIDENT  RAY  DIFFRACTED  RAY 


Figure  71.  Diffraction  at  Edge  of  Conducting  Wedge 

(Anderson,  1965:A-3) 


The  diffraction  coefficients  for  the  cylinder  shown  in 
Figure  72  can  now  be  derived  in  terms  of  the  incident  and 
diffracted  fields  angles  of  incidence  and  diffraction, 

and  WB,  respectively.  The  edges  contributing  to  the  RCS 
are  labelled  Sj,  S2,  and  S3.  The  fourth  corner  is  not 
illuminated  and  does  not  effect  RCS  if  the  cylinder  is  large 
compared  to  a  wavelength.  Conver'.ional  references  for  iDj 
and  are  shown  at  the  center  of  vhe  cylinder. 


I 

I 


Figure  72.  Diffraction  by  a  Right  Circular  Cylinder 

(Anderson,  1965:A-4) 


The  diffraction  coefficients  resulting  from  each  illuminated 
edge  in  Figure  72  are 


DlH  =  I  el  Tt/4S,in{>2/-3~  l  (Cos(2«/3)-Cos(tf8-^;)  ) 


-1 


( 2  irk ) 


-1 


T  (Cos  (  2  n/3  ) -Cos  (  2  (  K+tss+Vj ) /3  )  )  ]  (66a) 


D 


V  2  e**/4Sin(2n/3) 


2h  "  3 


(2*k) 


1/2^  [(C°s(2n/3)-Cos(^s-vi)) 


-1 


-vl 


T  (Cos  (  2n/3  ) -Cos  (  2  (  Fs+Fj ) /3  )  )  ]  (66b) 


P-U,  =  r--  —  4S" -^.g73-  [  (Cos(2-n/3)-Cos(^g-^j)  )_1 

3H  3  <2*k)1/2  8  1 


-1 


T  (Cos  (  2it/3  )  -Cos  (  2  (  «- Vg-Vj)  /  3  ) )  ']  (66c) 


The  diffraction  coefficients  given  in  eq.  (66)  are  for 
two-dimensional  right-angle  edges.  Because  a  cylinder  has 
curved  edges  it  is  three  dimensional.  The  curvature  is 
accounted  for  with  a  geometrical  spreading  factor  used  as  a 
weighting  function.  Letting  s  be  the  distance  from  the  edge 
to  the  receiver  along  the  scattered  ray  and  ^  be  the  radius 
curvature  of  the  of  the  diffracted  wave  front,  the 
weighting  function  can  be  written  as 


[s(l  +  Pj1®)]'172  (67) 

Because  s  >>  pj  ,  eq.(67)  can  be  written 

(  P!)1/2/s  as  Pj/s  — ►  0  (68) 
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is  related  to 


The  curvature  of  the  diffracted  wavefront,  j, 
the  angles  of  incidence  and  diffraction  and  the  cylinders 
curvature  by 


P1  = 


where 
Figure  73. 


Cos®  +  Cos6 

p  is  the  cylinder’s  radius  and  ©  and  6  are  defined 


(69) 


i  n 


Figure  73.  Angular  Relationships  for  Geometric  Factors 

(Anderson,  1965:A-7) 
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From  this,  the  geometrical  spreading  factors  for  each  edge 


are 

[(Pj1),/2]/s  =  t<P|2)1/2]/ s  r  { [a/Sinvj+SinVg)  l1/2)/s  (70a) 

[  (  Pi  3)  1/,21  /s  =  { [  a/  (  -S  invj  -  Sinv8)  ] 1//2} /s  (70b) 


The  minus  sign  occurs  in  eq.  (70b)  because  the  edge  at  S3 
curves  in  the  opposite  direction  when  compared  to  Sj  and  S2* 
The  phase  angle  associated  with  a  diffraction  point  is 
related  to  the  distances  the  incident  and  scattered  fields 
travel.  With  s  defined  as  the  distance  from  the  point  of 
diffraction  with  respect  to  a  fixed  reference  plane  and  b 
the  distance  the  incident  field  travels  relative  to  a  fixed 
reference  plane,  then  from  Figure  74,  s  =  rg  -  ds  ,  and 
b  =  rj  -  dj  ,  for  edge  Sj. 


Figure  74.  Phase  of  a  Diffracted  Ray 

(Anderson,  1965:A-8) 
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Tt  ran  be  derived  from  Figure  74  that 
dj  ^  =  aSinv'j  ♦  hCosttj 

and 

dR^  =  aSint>s  +hCos*>s 

Using  this  logic,  the  phase  terms  can  be  written 
ftik(S]+bj)  _  i  k  [  2i — a  ( Si  nt>j+R  invs ) -h  ( Cos«j+Cosvs  )  ] 

pik(s2+b2)  _  ftik[2i — a( Sin#>j+Sinvg)  +h ( Cos^+Cos^g)  ] 

pik(s^+b-^)  _  Rik  t  2r  +  a  (  Sin^j+S  i  nv>B) -h  ( Cos^j+Cosvs )  ] 


(71a) 


(71b) 


(72a) 


(72b) 


(  72o ) 


By  GTD,  RCS  as  a  function  of  the  angles  of  incidence 
and  diffraction  is  given  by 


o  {  Vj ,  VK ) 


4  nr^ 


2 

^total 

A 


(73) 


where  ^total  =  sum  aU  rays  diffracted  toward  the 

receiver 

Combining  eqs .  (66),  (70),  and  (72),  diffraction  at 
each  illuminated  edge  is  given  by 
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‘8  1 


2  Ae  />Sin(2n/3)  [«/ ( Sim»j4Sinv8)  ]  *'* 

** - 71^7?* - 5 - 

x  exp(ik[2r-a(Sini»j+Sin¥>8)-h(Cos¥>j+CoRtf8)  ] } 
x  [  (Cos(2n/3)-Cos(  V8-  Vj)  )_1 

*  ( Cos  (  2  n/3  )  -Cos  (  2  (  n+**s+i»j ) /3  )  )-1J 

2  Ae*  n/,<*Sin  (  2n/3  )  [a/ ( Sin^j+55 in vg )  ]  ^  ^ 


‘S2 


"3 


1/2 


(  2  nk ) 

x  exp{  i  k  [  2r-a(  SinVj+Sinp8)  +h  ( CosVj+Cosvs)  ] } 
x  [  (Cos{2n/3)-Cos(i>s-Vj)  f 1 

T  (Cos(2n/3)-Cos(2(^s+V.)/3)  f1] 


‘*3  =  *3 


2  Aein/4Sin(2«/3)  [a/ ( -Sin*rSin*>8)  ] 1/2 


1/2 


(74a) 


(74b) 


( Snk ) 

x  exp{ik[2r-a(SinVi+Sin^s)-h(CosVj+Cosvs)  ] ) 
x  [  (Cos(2n/3)-Cos(  Vg-Vj)  )_1 

T  ( Cos ( 2«/3 ) -Cos  (  2  (  n+ v8+t*i )  /  3  ) )~  *)  (74c) 


where  s  s  r  in  the  far  field. 

From  Figure  72  it  can  be  seen  vg  =  Vj  +  B  .  To  remain 
consistent  with  other  derivations  in  this  study  let  Vj  =  •  , 
Using  the  relationship  *>g  =  ©  +  B  ,  bistatic  RCS  can  be 
computed  as  a  function  of  the  angle  of  incidence  and  the 
bistatic  angle.  Eq.  (74)  is  now 


M 


»1 


2  Ae1 */4Sin(2n/3)  [«/ ( S i r«+Sin (•♦*)) 1/2 
(2nk)1/Z  r 

x  exp{ -i k [a ( Si n©+Sin ( ®+fi)  )  +h (Cos®+Cos  ( ®+g ) )  ] ) 
x  [ (Cos(2n/3)-Cos(2fi/3) )~J 

*  (Cos(2n/3)-Cos(2(n+2©+0)/3) )-1]  (75a) 


2  Ae1K/4Sin(2n/3)  [a/ ( S in©+Sin {©+*)] ^2 


M 


82  *  3  <2«k)1/2  7 

x  exp{-ik[a(Sin©+Sin(©+0) ) -h ( Cos©+Cos ( ©+0 ) )  ] ) 

x  [ (Cos(2n/3)-Cos(2fi/3) )_1 

*  ( Cos ( 2n/3 )  -Cos  (2(2©+j3)/3)  )-1] 


(75b) 


2  Aein/4Sin(2n/3)  t  a/ ( Sin©+Sin  (  ©+J5 )  ] 

Ms3  =  '3  F 

x  exp(ik[a(Sin®+Sin(®+£) ) -h ( Cos®+Cos ( ©+0)  )  J-in/2) 

x  [ (Cos(2*/3)-Cos(20/3) )-1 

*  (Cos(2n/3)-Cos(2(*-2©-g)/3) )_1]  (75c) 

Summing  eqs.(75a),  (75b),  and  (75c)  and  inserting  this  sum 
in  eq.(73)  gives  the  bistatic  RCS .  For  the  equivalent 
monostatic  RCS,  Vj  and  v>8  are  equal  to  ©  jB/2. 

Substituting  this  into  eq.(74)  gives 
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2  Ae1 "/4Sin(2*/3)  [a/ ( 2Sin ( •♦0/2 ) ) 1/2 

Msi  2  ^772  T 

x  exp{-i2k[aSin(©tf/2)+hCos(®+e/2)]} 
x  [ ( Cos ( 2  n/3 ) -Coa ( 20/3 ) )”4 

T  (Cos(2«/3)-Cos(2(*+2®+0)/3))'1]  (76a) 

2  Aein/4Sin(2*/3)  [a/ (  2Sin  ( ©tjS/2 )  ]  */2 
Mr2  =  7 

x  exp{-i2k[aSin(  e-ffl/2 )  -hCos(  ®+fl/2  )  ] } 

x  ( (Cos(2n/3)-Cos(2*/3) )-1 

T  ( Cos (2*/3) -Cost 2(2®+0)/3) )“*]  (76b) 

2  Ae1K/4Sin(2*/3)  [a/ ( 2Sin( e+0/2 ) ] 1/2 

Ws3  =  iw75  7 

x  exp{i2k[aSin(©+fl/2)-hCos(®+0/2) ]-i*/2) 
x  [ (Cos(2*/3)-Cos(2*/3) )_1 

*  (Cos(2x/3)-Cos(2(it-2e-fi)/3)  )~5]  (76c) 

Monostatic  RCS  is  then  confuted  by  summing  eqs .  (76a), 

(76b),  and  (76c)  in  eq.  (73). 

Tn  eqs.  (75c)  and  (76c)  the  phase  contains  the  term 
(-ix/2).  This  comes  from  (-Sin®j  -  Sin*>B)  in  the 
denominator  of  the  geometrical  spreading  factor. 
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non  non  onononono  nnnnn-nnn 


Computer  programs  written  in  FORTRAN  ??  to  calculate 
the  biatatic  and  equivalent  monostatic  RCS  are  listed  below 


C ******************<<<  BICYL . FOR  >>>*********************** 

THTS  PROGRAM  COMPUTES  THE  BISTATIC  RCS  OF  A 
RIGHT  CIRCULAR  CYLINDER,  HORIZONTAL  POLARIZATION. 

OUTPUT:  1  0*1.00 ( RCS  / 1  .  E4  )  ;  RCS  IN  SQUARE  CM 

*********:*  *********************************************** 

REAL  A , H , F , G , K , U , V , X , ALPH (451 ) .LAMBDA 
REAL  BET  A , PI , RC (461 ) ,S1 ,S2,S3,II ,Q,Q1 ,Q2,Q3 
COMPLEX  J , Y , PI , P2 , P3 
.1  =  (  0 .  ,  1  .  ) 

PT  =  3.141593 

DEFINE  VARIABLES 

A  =  RADIUS  OF  CYLINDER  LIN  CM 
H  s  HALF  CY  INDER’S  LENGTH  IN  CM 
LAMRDA  =  WAVELENGTH  IN  CM 
THETA  =  ANGLE  OF  INCIDENCE  (END-ON  =0*) 

BETA  =  BI STATIC  ANGLE 

A  s  7.50 
H  =  30.0 

OPEN  ( 1 ,STATUS='NFW’ , FILE= ’ B : SIGMA_BV . DAT ’ ) 

I  =  ANGLE  OF  INCIDENCE 

DO  40  I  =  0,90,15 

THETA  =  PI  * FLOAT ( I ) /l 8C . 

WRTTE ( 1 ,*)  ’THETA  =’  ,1 

N  =  BISTATIC  ANGLE 

DO  30  N  =  0,60 
WRITE ( * , * )  I , N 
BETA  =  FLOAT (N)*PI/180, 

LAMBDA  =  3.0 
K  =  2 . *PI/LAMBDA 

F  -  2.*SIN(2.*PI/3. )  /  (  3  . *SQRT ( 2 . *PI *K ) ) 

X  =  1 ./(COS(2.*PI/3. )-COS(2.*BETA/3. ) ) 

LI  =  180-N-I 
Q  =  COS ( 2 . *PI/3  .  ) 

Q1  s  COS ( 2 . * (PI+2 . *THETA+BETA ) / 3 . ) 

Q2  =  COS ( 2 . * ( 2 . * YHETA+BETA ) /3  .  ) 

Q3  =  COS ( 2 . * ( PI-2 . *THETA-BETA ) /3 . ) 

IF  (I.EQ.0)  THEN 
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c 

c 

c 

c 

c 


c 


20 

30 

C 

C 

C 


32 

35 

40 


RCS  s  PI*(A*t2*K*COS(BRTA/2. ) )**2 
ELSEIF  (I.EQ.90)  THEN 

RCS  s  '*K«A*H**2*COS(BETA/2. ) 


DTVTDE  BV  0  ERRORS  OCCUR  IN  THE  DIFFRACTION 
COEFFICIENT.  THE  FOLLOWING  INSTRUCTIONS  GIVE  THE 
RCS  FOR  THE  FREVTOUS  BTSTATIC  ANGLE  WHEN  THIS  OCCURS. 


El. SEIF  (Qt.EQ.Q)  THEN 
GO  TO  20 

ELSE  IF  (Q2.EQ.Q)  THEN 
GO  TO  20 

ELSE IF  (Q3.EQ.Q)  THEN 
GO  TO  20 

Rl.SETF  (I.EQ.L1)  THEN 
GO  TO  20 


ELSE 

G  =  SQRT( A/(SIN( ALPHA) +S IN (ALPHA+BETA) ) ) 

U  =  A*  ( SIN  (ALPHA)  4-SIN  (ALPHA+BFTA)  ) 

&  +H  * ( COS ( ALPHA ) +COS < ALPH A+BETA ) ) 

V  =  A»  (STN(  ALPHA) 4-SIN (AI.PHA  +  BETA)  ) 

A  -H* ( COS ( ALPHA ) +COS ( ALPHA+BETA ) ) 

'  '=  X+1 ./(Q-Ql ) 

S2=  X4- 1 . / ( Q-Q2 ) 

S3=  X+l ./(Q-Q3) 

Pl  =  F*GtS 1 *CEXP( J* ( PI/4 . -K*U ) ) 

P2s  F*G»S2*CEXP< J* (PT/4 . -K*V ) ) 

P3=  F*G*S3*CEXP(  J*(-l  ,#PT/4.4-K*V)  ) 

V  =  Pi +P2+P3 
RCS  =  4.*PT*CABS(Y)**2 
ENDTF 

RC ( N )  =  10.  *ALOG10(LCS/1  .E4  ) 

GO  TO  30 

M  =  N-l 

RC ( N )  s  RC ( M ) 

CONTINUE 

OUTPUT  THE  BISTATIC  ANGLE  AND  THE  CORRESPONDING  RCS 


DO  35  M  =  0,60 

WRITE (1,32)  M , RC ( M ) 
FORMAT ( I 3 , F 1 2 . 4 ) 
CONTINUE 
CONTINUE 
CLOSE (UNIT=1  ) 

END 


C 

C*********************************************************** 


c*t****tt**t*t*t ***<<<  MONOCYL . FOR  >>>*****#**»* ** ********* * 

C 
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onoo  odd  non  .onnonoonn  n  n  n  n  n  n  o  n 


THIS  PROGRAM  CALCULATES  THE  MONOSTATIC  RCS, 

HORIZONTAL  POLAARTZATION ,  CORRESPONDING  TO  A  GIVEN 
Bt STATIC  GEOMETRY  ACCORDING  TO  KELL’S  MONOSTATTC  C 
BISTATIC  RELATIONSHIP. 

OUTPUT:  10* LOG (RCS/ I .E4) ;  RCS  IN  SQUARE  CM 

9»*****»tt*M*»ttt****ttt*t**tt»t*t**t*t*«t**t*t*********t 

REAL  AIH,FIGIK,U,V,X,FREQ(451 ), LAMBDA 
REAL  BETA, PT,RC( 451 ) , SI ,S2 , S3 , IT ,Q,Q1 ,Q2 ,Q3 
COMPLEX  J , Y , PI , P2 , P3 
.1  =  ( 0  .  ,  1  .  ) 

PT  =  3.141593 

DEFINE  VARIABLES 

A  =  RADIUS  OF  CYLINDER  IN  CM 

H  =  HALF  CYLINDER’S  LENGTH  IN  CM 

LAMBDA  =  WAVELENGTH  IN  CM 

THETA  =  ANGLE  OF  INCIDENCE  (END-ON  s  0*) 

BETA  =  BISTATIC  ANGLE 

A  =  7.50 
H  =  30.0 

OPEN  ( 1 , STATUS* ’NEK’ ,FILE=  *B: SIGMA_FV.DAT’ ) 


I  =  ANGLE  OF  INCIDENCE 

DO  40  I  =  0,90,15 

WRITE (1 ,*)  ’ALPHA  =  ’  ,  I 

N  s  BISTATIC  ANGLE 

DO  30  N  =  0,60 

THETA  =  PI * ( FLOAT ( I ) +FLOAT (N)/2. ) / 1 80 . 

WRITE (*,*)  I , N 

BETA  =  FLOAT ( N ) *PI/180 . 

FREQUENCY  IS  ADJUSTED  ACCORDING  TO  KELL’S 
RELATIONSHIP. 

LAMBDA  =  3. 0*COS( BETA/2. ) 

K  =  2. *PI /LAMBDA 

F  =  2 . *SIN( 2 . *PI/3 . )/( 3 , *SQRT( 2 . *PI*K ) ) 
X  =  1 ,/(COS(2.*PI/3. )-l . ) 

THET  =  FLOAT(I)+FLOAT(N)/2. 

TERM  =  THET  ♦  FLOAT ( I ) 

Q  =  COS { 2 . *PI/3 . ) 

Q1  =  COS(2.*(PI*2.*THETA)/3. ) 

Q2  =  COS ( 2 . * ( 2 . *THETA ) /3 . ) 
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n  n 


Q3  s  COS( 2 . t ( PI-2 . 4THETA ) /3 . ) 

If  (THET.EQ.O. )  THEN 
RC.3  *  PII(AM2#K*COS(BETA/2. ) )M2 
ELSETF  (N.EQ.O)  THEN 
CO  TO  10 

10  ELSE IF  (THET.EQ.90)  THEN 

RCS  s  4*K»A*H**2*COS(BETA/2. ) 
ELSETF  (TERM. EQ. 210. )  THEN 
00  TO  SO 
ELSE 

G  s  SQRT (A/ (SIN ( THETA) ♦SINt THETA) ) ) 

U  s  A* (SIN (THETA) ♦SIN (THETA) ) 
t  +H M COS ( THETA )+COS( THETA ) ) 

V  =  A* ( SIN < THETA )+STN( THETA ) ) 

&  -H *( COS (THETA) ♦COS( THETA ) ) 

SU  X>1./(Q-Q1) 

S2=  X+l ./(Q-Q2) 

c1).  Y+1  /(Q.Q1) 

PI s  F*G*S1 tCEXPt J* { PI/4 . -K*U ) ) 

P2s  F*G*S2*CEXP( Jt ( PI/4 . -K*V) ) 

P3=  F#G#S3*CEXP(J4(-1 .*PI/4.+K*V) ) 

y  =  phps+ps 

RCS  =  4 .  *PI *CABS ( Y ) *  *  2 
ENDIF 

RC ( N )  =  1 0 . t AL0G1 0{RCS/1 . E4 ) 

30  CONTINUE 


OUTPUT  THE  BISTATIC  ANGLE  AND  CORRESPONDING  RCS 
DO  35  M  s  0,60 

WRITE( 1,32)  M,RC(M) 

32  FORMAT ( I 3 , F 1 2 . 4 ) 

35  CONTINUE 

40  CONTINUE 

CLOSE ( UNIT= 1 ) 

END 
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